das peças, ao contrário da transformação martensi- 
tica que, além de muito rápida logo que se atinge 
a temperatura M. no arrefecimento, não é uniforme 
em toda a espessura das peças visto a velocidade de 
arrefecimento ser sempre maior junto à superficie. 
Ora, para uma mesma temperatura, a densidade da 
austenite é maior do que a das estruturas em que se 
transforma, com ou sem difusão, pelo que há sempre 
aumento de volume nessas transformações, Tal au- 
mento origina tensões que, no caso da transformação 
martensítica, são obviamente muito mais perigosas 
devido à velocidade e heterogeneidade característica 
desta transformação, 

Os problemas de distorção e mesmo de fissuração 
das peças, especialmente quando de formas comple- 
xas, devido às variações de volume resultantes quer 
da dilatação ou contracção térmica, quer das diferen- 
cas de densidade das diversas estruturas, associadas 
à heterogeneidade das temperaturas através das pe- 
cas, são muitas vezes dificeis de resolver cabalmente, 
exigindo uma íntima colaboração entre os projectis- 
tas e a secção de tratamentos térmicos, 

Também a verificação de eventuais fissuras, su- 
perficiais ou interiores, bem como da profundidade 
de endurecimento obtida, exige todo um conjunto de 
ensaios não destrutivos, desde os líquidos penetrantes 
aos ultra-sons e correntes de Foucault, cuja interpre- 
tação reveste por vezes apreciável dificuldade. 

Ainda outres problemas também de grande im- 
portância são os ocasionados pelas operações de recti- 
ficação final das peças, pelos perigos que acarretam 
de aberturas de fendas superficiais devido às tensões 
residuais ou de amaciamento local por aquecimento 
exagerado, bem como a protecção das peças contra 
a descarburação superficial durante a austenitização 
(que obriga muitas vezes à utilização de banhos de 
sais ou de atmosferas controladas), e o controlo 
rigoroso das temperaturas de tratamento. 


As estruturas resultantes tanto da austêmpera 
como do revenido da martensite dependem da tempe- 
ratura a que foram obtidas, pelo que é possível con- 
seguir uma larga gama de propriedades para os 
acos tratados. No entanto, a previsão quer da estru- 
tura que resultará de dado ciclo térmico quer das 
propriedades que dessa estrutura resultarão, embora 
possa fazer-se com uma certa aproximação, mercê 
da cada vez maior quantidade (e qualidade) de in- 
formação que vai sendo acumulada, não evita que se 
mantenha ainda uma larga margem de incerteza, 
obrigando a investigação caso a caso para uma opti- 
mização do tratamento, 

Tem sido praticamente apenas no estrangeiro 
que se têm desenvolvido os conhecimentos básicos 
nesta área. À par de estudos cada vez mais porme- 
norizados sobre as diversas microestruturas e a sua 
génese, bem como sobre as propriedades que delas 
resultam (através da utilização de microscopia de luz 
visível e electrónica de transmissão e de varrimento, 
de ensaios de difracção, dilatométricos e mecânicos, 
de análise química por micro-sonda electrónica, etc.), 
têm-se desenvolvido novos métodos de tratamento fa- 
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zendo uso dos progressos da ciência e da técnica 
noutros domínios, cs quais permitem uma posição de 
vanguarda dos paises neles desenvolvidos. 

Assim, por exemplo, dado que em muitos casos 
a propriedade mais procurada é uma boa resistência 
ao desgaste, promovem-se muitas vezes modificações 
estruturais apenas numa fina película superficial das 
peças. Estas modificações podem conseguir-se quer 
por simples têmpera superficial quer por alteração 
da composição química dessa película. No primeiro 
caso austenitiza-se apenas a película a endurecer, 
por meio de um aquecimento muito rápido através 
de correntes eléctricas induzidas de alta frequência 
aproveitando o «efeito pelicular>», ou por incidência 
na superficie de raios «laser», de chamas como a 
oxi-acetilênica, e ainda por bombardeamento de elec- 
trões sob vácuo, seguindo-se um arrefecimento tam- 
bém rápido que, conforme os casos, pode ser por 
transmissão do calor para um meio exterior às peças 
ou para o próprio interior, ainda frio, das peças. 
No segundo caso, em que se incluem a cementação, 
a nitruração, a boretação, etc., usam-se também mé- 
todos cada vez mais sofisticados por forma a serem 
conseguidas as estruturas de muito alta dureza (com 
a boretação conseguem-se durezas da ordem de 2 000 
unidades Vickers) fazendo o tratamento apenas nas 
zonas desejadas, evitando a oxidação das superfícies, 
distorcendo o mínimo possível as peças, e tornando 
os processos cada vez mais seguros para os opera- 
dores. 

Têm igualmente sido desenvolvidos trabalhos de 
investigação no sentido de conferir aos aços a melhor 
estrutura para uma maquinabilidade ou uma capaci- 
dade de enformação óptimas, o que se pode conseguir 
também muitas vezes por tratamentos térmicos ade- 
quados. 

A finalizar este capítulo, no qual se procurou 
chamar a atenção para os principais tipos de proble- 
mas que surgem nos tratamentos térmicos dos aços, 
faz-se notar qce continuam a desenvolver-se estudos 
quer para uma melhor elaboração dos aços já exis- 
tentes (principalmente no sentido de lhes conferir a 
maicr homogeneidade em toda a espessura dos per- 
fis), quer para a criação de novos aços (é o caso, 
por exemplo, da criação ainda recente dos aços «ma- 
raging» que são ligas de Ferro e Níquel a que se 
juntam outros elementos tais como o Cobalto, o Mo- 
libdénio, o Alumínio, o Titânio e o Vanádio, sendo 
o teor em Carbono o mais baixo possível; nestes aços 
a transformação martensítica consegue-se mesmo em 
arrefecimento lento, promovendo-se seguidamente o 
endurecimento através da precipitação muito fina de 
compostos intermetálicos no seio da martensite por 
envelhecimento (aging) desta, e daí o nome dos aços). 


2-— OS TRATAMENTOS TÉRMICOS DOS AÇOS EM 
PORTUGAL 


Portugal importa anualmente várias dezenas, se 
não centenas, de milhares de contos de aços destina- 
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dos ao fabrico de peças a serem tratadas termica- 
mente. Estes aços, pela sua cuidada elaboração e com- 
posição, são necessariamente caros, e tal facto, asso- 
ciado às grandes despesas em equipamento e mão- 
-de-obra que o fabrico das diversas peças acarreta, 
justifica plenamente que se encare o problema do 
tratamento térmico com o maior cuidado, Infelizmen- 
te, muito poucas são no nosso País as empresas que 
dispõem de instalações convenientemente apetrecha- 
das para o tratamento térmico dos aços que empre- 
gam e de pessoal devidamente preparado para esse 
efeito. 

De um modo geral, os tratamentos são feitos 
com base apenas nas indicações de carácter gené- 
rico que os próprios fornecedores dos aços facultam, 
sendo raríssimos os casos em que se efectuam estu- 
dos com vista a optimizar os tratamentos. Não admira 
portanto que muitas peças tenham de ser rejeitadas 
logo após o tratamento térmico ou venham a ter du- 
ração muito curta em serviço, resultando daí grandes 
perdas para o Pais. Por outro lado não existe também 
ainda uma coordenação de esforços devidamente es- 
truturada entre fornecedores e utilizadores dos aços 
de tratamento térmico por um lado e organismos de 
investigação por outro, capaz de permitir um am- 
biente de troca permanente de informações e de 
incentivação de resolução conjunta dos problemas. 


Alguma coisa tem sido feita, no entanto, que 
poderá servir de base para Portugal se situar rapi- 
damente ao nível de outros países nos quais associa- 
ções especializadas, apoiadas em laboratórios devi- 
damente equipados, incentivam continuamente o pro- 
gresso nesta área promovendo congressos, Cursos, 
seminários, publicações, periódicas ou não, etc. Uma 
ou outra publicação, de carácter genérico na maioria 
dos casos, tem sido editada quer por entidades públi- 
cas quer privadas, e uma ou outra actividade forma- 
tiva tem sido também realizada, As actividades mais 
recentes consistiram na realização de um ciclo de 
seminários sobre tratamentos térmicos dos aços pro- 
movido pelo Grupo de Estudos da Tecnologia dos 
Aços de Construção Mecânica da Faculdade de Ciên- 
cias e Tecnologia da Universidade de Coimbra, que 
teve a colaboração do Laboratório Nacional de Enge- 
nharia Cívil e do Centro de Mecânica e de Materiais 
das Universidades de Lisboa — CEMUL, num curto 
curso de tratamentos térmicos dos aços da iniciativa 
do Instituto de Soldadura, e num Curso de Especia- 
lização em Ciência dos [Materiais organizado pela 
Universidade Nova de Lisboa. O CEMUL, recente- 
mente criado, incluiu também nas suas linhas de 
acção o estudo dos tratamentos térmicos dos aços, 
tendo a decorrer um projecto de selecção de acos- 
“ferramenta com vista à sua aplicação à indústria 
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nacional em colaboração com o Laboratório de Me- 
talurgia do Instituto Nacional de Investigação Indus- 
trial. 

Há por conseguinte um certo número de pessoas 
interessadas em promover o desenvolvimento neste 
campo, sendo do maior interesse associá-las desde 
já para realizarem um trabalho conjunto. Este grupo 
seria de esperar que viesse a aumentar rapidamente 
pela adesão das pessoas que nas diversas empresas 
se preocupam com estes assuntos, bem como de téc- 
nicos das firmas fornecedoras quer dos aços quer do 
equipamento de tratamento térmico, Dado que as ope- 
rações de soldadura dos aços implicam ciclos térmi- 
cos, é de esperar que as pessoas interessadas na sol- 
dadura também se venham a juntar a este grupo. 

O autor da presente comunicação sugere que to- 
das as pessoas interessadas sejam postas em con- 
tacto, de início, com base na Ordem dos Engenheiros, 
até se criar uma associação técnica autónoma a 
exemplo de outras nações, e sugere também que entre 
as tarefas primordiais a efectivar se inclua a incen- 
tivação da criação pelo Governo de uma oficina 
nacional de tratamentos térmicos dos aços, de prefe- 
rência integrada no previsto Laboratório Nacional de 
Engenharia e Tecnologia Industrial, do Ministério da 
Indústria e Tecnologia, bem como da dotação do 
equipamento de investigação indispensável aos orga- 
nismos que em Portugal têm por missão a investiga- 
cão no âmbito da Ciência dos Materiais, 

Tal oficina nacional teria como missão: suprir 
a falta de instalações adequadas ao tratamento tér- 
mico dos acos naquelas empresas públicas ou pri- 
vadas que delas careçam, e em que muitas vezes não 
é rentável terem-nas, sendo preferível o recurso a 
meios externos; apoiar as empresas que dispõem de 
instalações de tratamentos térmicos na resolução dos 
seus problemas, quer efectuando estudos de optimiza- 
cão de determinados tratamentos, quer colaborando 
na manutenção de níveis adequados de trabalho atra- 
vés de cursos, estágios, etc. destinados ao pessoal 
dessas empresas, de vigilância do bom funcionamento 
do equipamento de medida e controlo, etc.; transmi- 
tir a organismos universitários de investigação aque- 
les problemas que transcendam os seus meios, dando 
a esses organismos o apoio possível na investigação 
a efectuar; criar e manter actualizados ficheiros 
com os quantitativos e custos de todos os aços im- 
portados, com base nos Serviços Aduaneiros Nacio- 
nais, por forma a racionalizarem-se as importações 
em presença das utilizações dadas no Pais a esses 
aços, podendo tal ficheiro, conjuntamente com a ra- 
cionalização referida, vir a permitir uma futura to- 
mada de decisão sobre a eventual criação de uma 
acearia de aços para tratamento térmico em Portugal, 
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A aglomeração de concentrados finos de minérios 


RESUMO 


Nos modernos processos da metalurgia extracti- 
va exigem-se cargas devidamente apropriadas de 
modo a conseguirem-se produtividades tão elevadas 
quanto possível e custos energéticos tão reduzidos 
quanto possível. 

Neste trabalho indicam-se os processos de aglo- 
meração mais difundidos mundialmente e que tem a 
sua maior aplicação na indústria siderúrgica: sinte- 
rização e peletização. 

Trata-se depois de inventariar o que já existe 
em Portugal e quais as perspectivas futuras tendo 
em atenção os recursos nacionais em minérios de 
ferro. 


1 — INTRODUÇÃO 


A necessidade de proceder ao enriquecimento 
dos minérios naturais, por razões de ordem econó- 
mica levou à necessidade de instituir uma operação 
pré-metalúrgica que se designa por aglomeração. Com 
efeito, para se enriquecer há que proceder a fragmen- 
tações e por isso mesmo os concentrados serão pro- 
dutos de pequena granulometria que para serem car- 
regados nos fornos de redução terão de primeiramente 
ser aglomerados, 

Embora normalmente se liguem as operações de 
aglomeração de concentrados à indústria siderúrgica 
uma vez que é aqui que as instalações atingem dimen- 
sões espectaculares e a técnica se encontra muito 
avançada, a aglomeração está generalizada e aplica-se 
à preparação das cargas para redução utilizando 
aparelhos com dimensões muito diversas. Assim, en- 
contram-se instalações de aglomeração na apropria- 
ção de concentrados de sulfuretos (galenas, blendas, 
etc.) ao mesmo tempo que se promove a ustulação à 
morte, instalações essas por vezes apenas para algu- 
mas centenas de toneladas por ano, até instalações 
de muito grande dimensão para a alimentação de al- 
tos fornos que podem produzir algumas dezenas de 
milhar de toneladas por dia, 
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ABSTRACT 


In extractive metallurgy modern processes for 
preparing the charges are required in order to reach 
as high as possible productivity and as low as possible 
energetic costs. 

In this work are indicated the more worldlhy 
widspread agglomeration systems and those that have 
their main application in the siderurgical industry: 
sintering and pelletizing. 

After that an inventory of what already exists 
in Portugal is draum up together with its future 
prospects, paying attention to the national resources 
in iron ore. 


E na aglomeração de concentrados de minérios 
de ferro que se encontram as instalações gigantes 
embora noutros domínios se assista a um certo pro- 
gresso nomeadamente na aglomeração de minérios de 
manganês. 

Dentro da aglomeração de minérios de ferro a 
que dedicaremos o fundamental deste trabalho exis- 
tem dois processos hoje generalizados no mundo e 
que são a sinterização e a peletização. Embora a 
peletização seja o processo industrial mais recente 
tem tido um desenvolvimento muito rápido mercê de 
por um lado haver a necessidade de enriquecer cada 
vez mais os produtos a carregar nos altos fornos a 
fim de diminuir o consumo específico de coque e por 
outro irem desaparecendo os minérios de dessimina- 
ção grosseira que obrigam, para que a libertação dos 
óxidos de ferro tenha lugar, a moagens extremamente 
finas. A aglomeração por sinterização, mais antiga 
no domínio industrial, tem sido objecto de grandes 
aperfeiçoamentos mercê de uma investigação profun- 
da pelo que se manterá, enquanto dispuser de miné- 
rios com dimensões apropriadas, com produções supe- 
riores. Isto se deve ao facto de não só haver um 
grande número de instalações em produção mas tam- 
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bém porque os <«altofornistas» consideram que os 
e«sintera» se mantêm na passagem através do alto 
forno mais uniformemente distribuídos, não se auto- 
-segregando devido à sua forma extremamente irre- 
gular. 

Será, no entanto, inevitável a substituição dos 
«sinters» pelas peletes à medida da exaustão dos mi- 
néios ricos obrigando à utilização de minérios pobres 
de dessiminação fina. 


Porém, a peletização poderá vir a substituir a 
n sinterização mais cedo do que se pensa, Com efeito, 
um pouco por toda a parte no mundo industrializado 
sº investiga a possibilidade de passar da siderurgia 
clássica (alto-forno, acearia a oxigénio) para uma 
siderurgia nova onde o alto-forno por razões eco- 
nómicas seja eliminado ficando em seu lugar apare- 
lhos de redução no estado sólido, Os aparelhos que 
melhor se adaptam a este trabalho de pré-redução 
continuam a ser fornos de cuba, reduzidos à cuba 
propriamente dita e onde a carga é constituída exclu- 
sivamente por peletes e a redução se efectua com 
recurso a gases redutores. Quem mais terá a ganhar 
com o desenvolvimento desta nova siderurgia são os 
países em vias de desenvolvimento que possuindo com 
frequência minérios, gases naturais, produtos deriva- 
dos do petróleo ou carvões não coquizáveis, que não 
consomem, poderão instalar mini-aceraias eléctricas 
para a produção do aço de que necessitam alimen- 
tando-as com peletes pré-reduzidas. 

Estas poderão ainda substituir sucatas de aço 
com vantagens em muitos casos, principalmente por- 
que, com facilidade, conduzirão a um produto de 
grande homogeneidade ao longo do tempo e sem o 
perigo da introdução de elementos por vezes nocivos 
que as sucatas com frequência contém. 


O trabalho que iremos desenvolver em seguida 
tem como finalidad: dar a conhecer técnicas que se 
generalizaram no mundo siderúrgico e uma experiên- 
ci2 pessoal adquirida ao longo do tempo quer em 
trabalho industrial, trabalho laboratorial ou visitas 
a instalações de pequena e grande capacidade. 


2 — SINTERIZAÇÃO 


Desde a instalação Greenwalt de Vila Cova com 
uma capacidade de cerca de 120 tons/dia de «sinter> 
até à instalação Dwight Lloyd de Dunquerque para 
uma produção de 16 000 tons/dia conhecem-se cente- 
nas de instalações com produções muito diversas. 
Todas se situam ma proximidade do aparelho de re- 
dução a que os «sinters» se destinam uma vez que 
estes não podem ser transportados sem sofrerem 
degradações consideráveis. Com efeito, trata-se de 
um material muito duro parcialmente vitrificado que 
se degrada por abrasão, reduzindo-se a partículas de 
pequenas dimensões. 
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Embora fixos, semi-fixos ou móveis os aparelhos 
cujas hases são constituídas por grades de barras, 
são carregados com uma carga de minério, castinas 
e combustíveis de altura previamente definida obri- 
gando-se a chama, após ignição com queimadores 
apropriados, a atravessar toda a carga, alimentando- 
-se a combustão com ar aspirado inferiormente. Ape- 
sar de a quantidade de combustível ser pequena (4 a 
E % em peso da carga mineral) as temperaturas atin- 
gem 1400 e mesmo 1600" €C, o que é suficiente para 
levar a carga a uma fusão incipiente o que provoca, 
pelo menos, a «colagem»» das partículas entre si. 

A operação de sinterização é por isso uma ope- 
ração de combustão de um combustível húmido com 
elevado teor em cinzas, usando vento como combu- 
rente, produzindo elevadas temperaturas, com obten- 
ção de um produto poroso semelhante a uma escória, 
denominado sínter. 

Para além da combustão haverá uma transferên- 
cia de calor que se transmite da carga para o ar e 
deste de novo para a carga de modo a que, se tudo 
estiver bem regulado, a temperatura de ignição do 
combustível será atingida quando estiverem perfeita- 
mente estabelecidas as condições de comburência. 

Deste modo, devido à transmissão da combustão 
de cima para baixo, através da camada a aglomerar, 
e devido também à diminuição das perdas de calor 
para o exterior através das paredes do aparelho e 
através da zona superior já aglomerada, a tempera- 
tura da zona de combustão vai subindo à medida que 
o tempo decorre e só não atinge o valor máximo 
ao chegar à grelha de suporte porque esta é protegida 
por uma camada de aglomerado com 3 a Sem de 
altura (produzido numa operação anterior, de gra- 
nulometria conveniente) e com uma capacidade calo- 
rífica suficiente para que a temperatura da grelha 
não ultrapasse os 500º C e se não deteriore rapida- 
mente, 

A investigação laboratorial dos fenómenos em 
jogo levou à conclusão de que haveria que aspirar 
um dado volume de ar, para se conseguir a aglome- 
ração de uma dada carga, o qual era determinado não 
pelas necessidades da combustão mas sim pelas de 
transferência de calor. Assim se pode, em primeira 
aproximação, considerar a existência de um volume de 
ar específico para a aglomeração cujo peso é sensivel- 
mente igual ao peso da carga a aglomerar. Tal facto 
pressupõe que os calores específicos da carga e do 
fluído aspirado sejam iguais o que não está longe da 
realidade. 

Para que se possa promover a aspiração do ar 
necessário à operação num dado tempo, compatível 
com uma boa aglomeração, torna-se necessário que 
a carga tenha elevada permeabilidade ao ar e que 
seja por este atravessada quando se estabelece uma 
depressão inferiormente à grelha usando para o efeito 
um ventilador, O tempo para uma operação de aglo- 
meração deverá ser determinado à «priori» e a partir 
do seu conhecimento poderá dimensionar-se a Insta- 
lação a construir. 

Das razões apontadas ressalta o facto de haver 
necessidade imperiosa de cuidar da apropriação da 
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* Equipamentos para cimenteiras 
x Equipamentos para siderurgias 
* Equipamentos para centrais de 
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pecificações API 
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carga a lançar sobre o aparelho de aglomeração, 
apropriação essa que terá de ter em atenção a na- 
tureza e a composição granulométrica das substân- 
cias minerais (minérios, castinas, retornos), a compo- 
sição pgranulométrica e a quantidade dos finos de 
coque e ainda o teor em humidade e a altura da ca- 
mada, 

A modificação de qualquer destes factores impli- 
cará a modificação dos restantes a fim de procurar 
manter as características físicas da carga e conse- 
quentemente o tempo de aglomeração. Por isso mes- 
mo é que a modificação da natureza e/ou granulo- 
metria do minério que é o componente principal da 
carga poderá ter uma influência enorme no decurso 
da operação, traduzida pelo tempo que a combustão 
leva a transmitir-se da zona superior à zona inferior, 
A diminuição da permeabilidade poderá obrigar a 
modificar a altura da camada e consequentemente da 
capacidade de produção da instalação. 

Dado que, tal como dissemos mais acima, a ten- 
dência é para a produção de concentrados cada vez 
mais finos e como se torna bastante oneroso modi- 
ficar as características das instalações, multiplicam-se 
os estudos laboratoriais conducentes à mlhoria da 
permeabilidade das cargas a aglomerar. No fundo o 
que se procura é modificar um factor até certo ponto 
subjectivo e que os franceses intitulam como «aptidão 
à aglomeração». No fundo a «aptidão à aglomeração» 
dependerá da capacidade de as partículas mais finas 
formarem umas com as outras e com as de maior 
calibre fragmentos de maiores dimensões dando lugar 
a um empilhamento em que as partículas finas não 
colmatem os espaços existentes entre as partículas 
grosseiras. Tudo se deverá passar como se houvesse 
uma classificação mais cerrada das partículas cons- 
tituitivas da carga. A melhoria da «aptidão à aglome- 
ração» passará em primeiro lugar pelo doseamento 
correcto da humidade, por uma operação de mistura 
em aparelhos devidamente adaptados à obtenção dos 
efeitos enunciados e em última análise pela adição 
de um ligante à semelhança do que se faz em pele- 
tização. 

Nas instalações siderúrgicas integradas usa-se 
com frequência como ligante o sulfato ferroso re- 
circulado a partir de instalações de decapagem da 
chapa e também o «schlames> de lavagem dos gases 
dos convertidores de acearia. Os principais inconve- 
nientes advêm da adição de S e compostos alcalinos. 

No caso de estes produtos não poderem ser usa- 
dos procede-se à substituição total ou parcial da cas- 
tina por cal. Autores brasileiros conseguiram deste 
modo aumentar a capacidade de uma dada instalação 
em cerca de 40 GG. Nas instalações laboratoriais da 
Faculdade de Engenharia foi possível por utilização 
de uma técnica semelhante melhorar a permeabi- 
lidade de uma camada de altura igual a Sem, 
obtida com minérios magnetíticos do Marão muito 
finos (70 % — 200 mesh) e de baixa aptidão à aglo- 
meração, diminuindo assim o tempo da operacão. 
Na transposição para a prática industrial foi possi- 
vel aumentar a capacidade de instalação de aglo- 
meração existente em cerca de 2064. No entanto 
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para a operação prévia de mistura não foi possivel 
montarem-se os meios necessários e mais convenien- 
tes e daí a ocorrência de uma certa irregularidade 
de marcha, Na prática a melhoria considerável só 
poderia vir a obter-se se se conseguisse uma verda- 
deira micropeletização, num misturador devidamente 
adaptado, do produto a aglomerar. 

A composição da carga deverá ser estabelecida 
tendo em ateação a obtenção de um leito de fusão 
para o alto-forno. Assim, da carga deverão fazer 
parte SiO,, CaO, AL,O, e MgO convenientes para o 
estabelecimento de uma escória qualitativa e quan- 
titativamente adaptada à marcha do alto-forno. Os 

Cao + MgO 


SiO, + ALO, 
deverão ser aproximadamente 1,2 e 0,9 respectiva- 
mente, 

A temperatura a atingir-se durante a operação, 
deverá ser suficiente para a obtenção de um aglome- 
rado resistento às acções de degradação tanto no 
exterior como no interior do alto forno mas, não tão 
elevada que provoque uma fusão total com formação 
de um produto compacto de baixa superficie especi- 
fica e fundamentalmente pouco redutível. A redutibi- 
lidade dependerá fundamentalmente dos constituintes 
que resultarem da dissociação e combinação dos com- 
ponentes da carga, Sob este ponto de vista haverá 
que ter em conta os constituintes que se formam 
não só por reacções no estado sólido como também 
aqueles que resultam da solidificação após fusão mais 
ou menos incipiente dos componentes. 

A diminuição do potencial químico de oxigénio 
determinará dificuldades de redução o que se tradu- 
zirá por aumento do consumo de coque para a redu- 
ção no alto forno. Por isso mesmo a temperatura 
deverá situar-se entre limites bem definidos sem o 
que ou o sinter se degradará facilmente ou se tor- 
nará pouco redutível. 


índices de basicidade CaO/SiO, ep 


3 — PELETIZAÇÃO 


Muito diferente da sinterização, a peletização foi 
de início uma operação que procurava fabricar a par- 
tir de concentrados ricos e de calibre muito fino um 
produto de dimensões suficientes para carregamento 
directo nos fornos de redução. Na prática os produtos 
obtidos poderiam ser utilizados como se de minérios 
siliciosos e ricos se tratasse. Portanto, e diferente- 
mente da sinterização, as peletes não poderiam por 
si só constituir um leito de fusão pois como possuíam 
um índice de basicidade muitíssimo baixo dariam 
lugar à produção de escórias muito viscosas e sem 
a capacidade depuradora que normalmente se lhes 
exige, 

Actualmente ainda a maior parte das peletes 
produzidas são siliciosas tornando-se por isso necessá- 
rio constituir leitos de fusão que as têm por base 
adicionando castinas ou sinters básicos. No entanto, 
em certas circunstâncias vão sendo produzidas pele- 
tes autofundentes que podem substituir total ou par- 
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cialmente os «sinters» e peletes pré-reduzidas com que 
se pretende enriquecer os leitos de fusão dos altos 
fornos, ou substituir total ou parcialmente as sucatas 
seja nos convertidores, como agente arrefecedor, seja 
nos fornos eléctricos, para a produção directa de aço. 

O fabrico de peletes seja de que tipo for decorre 
em várias etapes bem distintas: preparação da carga 
por mistura dos constituintes (minérios e ligantes), 
fabrico das peletes em aparelhos de peletização pro- 
priamente ditos de que se destacam os discos e os 
tambores por serem até agora os mais usados, € co- 
zedura das peletes com a finalidade de promover o 
endurecimento necessário ao seu transporte no ex- 
terior e no interior dos fornos. Para efectuar a coze- 
dura os aparelhos mais usados são os fornos de cuba, 
as cadeias tipo Dwight-Lloyd, semelhantes aos apare- 
lhos de sinterização e os fornos «grate-kiln», mistos 
de cadeia e forno rotativo. 

Não vamos fazer o estudo de cada uma das eta- 
pes de fabrico uma vez que se trata de operações 
bem conhecidas dos engenheiros que se interessam 
por estes problemas. Interessa no entanto dizer que 
a finalidade é aglomerar partículas finas (80 a 90 % 
— 44 micra) de modo a obter bolas com um diâmetro 
médio de 10/12 mm, quando destinadas à carga do 
alto forno. Quando destinadas a acearia de conversão 
o diâmetro poderá crescer até cerca de 30mm uma 
vez que se destinam a substituir fragmentos de miné- 
rio rochoso, rico e denso. O fabrico das peletes resulta 
de um processo muito parecido com o da formação 
de bolas de neve. Assim, gotoletas de água caídas 
nos peletizadores dão lugar a núcleos que crescem 
à custa de outros ou de particulas secas sobre que 
rolam segundo mecanismos onde a tensão superficial 
da água desempenha papel relevante, 

A cozedura tem por finalidade não a sinterização, 
isto é, a fusão das particulas e sua coalescência, mas 
sim uma recristalização, com ou sem oxidação, e for- 
mação de compostos no estado sólido. Daí que os 
constituintes observados nas peletes cozidas sejam 
muito diferentes dos observados nos +sinters», Assim, 
nas peletes são criadas pontes entre os grãos de óxido 
de ferro as quais são facilitadas se entretanto houver 
oxidação (caso em que a magnetite passe a hematite) 
ou são formados compostos entre os componentes das 
gangas ou dos ligantes e os próprios óxidos de ferro. 
De qualquer maneira os constituintes que envolvem 
elementos de ganga tem que formar uma rede con- 
tíinua muito estável até à fusão. Com efeito, a resis- 
tência das peletes à compressão terá que ser encarada 
não só durante a sua movimentacão no exterior dos 
fornos de redução mas principalmente no decurso da 
redução no seu interior. Na realidade a redução dos 
óxidos de ferro é realizada em grande parte no estado 
sólido pelos gases redutores CO e H, dando lugar à 
formação de glóbulos de ferro muitas vezes isolados 
no seio das peletes. A rotura das pontes entre par- 
tículas de óxido de ferro promoveria a diminuição 
da resistência das peletes aos esforços de compressão 
e abrasão (atrito) a que estão sujeitas quando a tem- 
peratura se situa entre os 1000 e os 1200º € em meto 
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redutor. Daí que pudesse ocorrer o seu esmagamento, 
redução a pó fino e seu transporte pelos gases. Nessas 
zonas do forno desempenha papel fundamental a liga- 
cão escória que só desaparece quando da fusão, Por 
isso, o teor em ferro das peletes terá um limite supe- 
rior a não ultrapassar e que deverá ser determinado 
em ensaios laboratoriais de redutibilidade. 

As peletes autofundentes tém em geral um teor 
era ferro mais baixo que as peletes ácidas fabricadas 
com o mesmo minério pois, haverá que adicionar uma 
quantidade de ligantes suficiente à base de Ola para 
o estabelecimento de um índice de basicidade conve- 
niente sempre muito mais elevado que o das peletes 
ácidas, 

A cozedura destas peletes é mais difícil uma vez 
que terá de decorrer emtre limites de temperatura 
mais apertados, determinados pela maior facilidade 
de fusão dos constituintes formados a partir da SiO, 
e Ca0O. 

As peletes pré-reduzidas começaram por ser na 
maior parte dos casos peletes ácidas que são em se- 
guida reduzidas em aparelhos onde circulam gases 
redutores a temperaturas convenientes mas que nor- 
malmente não ultrapassam os 1200" C. No entanto, 
outrós processos são aplicáveis com recurso a com- 
bustiveis sólidos e líquidos, Os combustiveis sólidos 
podem ou ser misturados com as peletes e aquecidas 
em conjunto a fim de se produzir a atmosfera con- 
veniente para, à temperatura de trabalho, se realizar 
a redução, ou serem mesmo incorporados na mistura, 
previamente à peletização, ocorrendo a pré-redução 
quando da cozedura, Experiências laboratoriais per- 
mitem antever um progresso notável nestas técnicas. 

Por fim não queremos deixar de referir que na 
evolução técnica da peletizacão surgiu recentemente 
o processo Grangcold. Utilizado na aglomeração de 
minérios de magmetite de mina sueca de Grangsberg 
no processo Grangecold hã apenas a considerar as fa- 
ses de mistura e de peletização. Com efeito, o pro- 
cesso surgiu com a finalidade de evitar a fase mais 
cara que é justamente a cozedura. Neste processo o 
ligante pode ser cimento Portland ou cal hidráulica 
em vez de bentonite, No entanto, a quantidade de 
ligante que no caso de utilização da bentonite atin- 
girá no máximo 1 % aqui pode atingir 10 %. Por ou- 
tro lado as peletes terão de ganhar presa, tal como 
acontece nas obras de betão, e por isso deverão ser 
deixadas em parque durante pelo menos quatro dias. 
Só após esse tempo de repouso, durante o qual se dá 
o endurecimento com um aumento de resistência à 
compressão notável, elas serão susceptíveis de mani- 
pulação, 

Este processo tem do ponto de vista econômico 
enorme interesse uma vez que os investimentos para 
a produção deste tipo de peletes são muito inferiores 
aos das instalações normais. No entanto, parece que 
os inconvenientes são bastantes para não permitir, 
pelo ménos a curto prazo, a sua generalização. De 
acordo com algumas informações que conseguimos co- 
lher o inchaço que sofrem estas peletes em curso de 
redução é muitissimo elevado o que pode provocar 
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acidentes de marcha semelhantes aos causados pela 
subida do teor em elementos alcalinos na carga. 

Daí que a carga dos altos fornos não tenha ainda 
podido contar com mais de 15 % de peletes grangcold, 

É curioso notar que peletes deste tipo foram 
talvez usadas industrialmente pela primeira vez em 
Portugal entre 1960 e 1967. Na realidade após uma 
fase inicial em que o forno Tysland Hole de Vila Cova 
do Marão usou um leito de fusão frio constituído por 
sinter, procedeu-se à sua transformação que no essen- 
cial consistiu na montagem de chaminés para carga, 
pré-aquecimento e pré-redução. Simultaneamente e 
porque se pretendia obter uma gusa de mais baixo 
teor em fósforo foi preciso purificar o minério moen- 
do-o mais fino e reconcentrando-o. Como a sinteriza- 
cão se não adaptava a este novo tipo de concentrado 
foi necessário montar uma instalação de peletização. 
A mistura para a obtenção das peletes era constituída 
por minério, carvão finamente moído e cimento que, 
na proporção de cerca de 10 %, actuava como ligante. 
Estas peletes eram também armazenadas durante 
4 ou mais dias a fim de endurecerem e virem a ser 
em seguida pré-aquecidas e pré-reduzidas ao cami- 
nharem para o forno. É curioso notar que este ais- 
tema raramente funcionou correctamente devido ao 
encravamento das chaminés de carga. Embora se não 
tivesse determinado com rigor as razões do encrava- 
mento hoje pedemos suspeitar ter sido o inchaço das 
peletes que o provocava. 

Actualmente produzem-se cm Vila Cova peletes 
grangeold usando como ligante cimento na proporção 
de 8 %. Estas peletes após repouso de 5 a 7 dias pode- 
rão ser manipuladas e carregadas no forno eléctrico 
sem problemas de degradação notável, 


4 — AGLOMERAÇÃO EM PORTUGAL 
As instalações de que se dispõe actualmente são: 


a) Sinterização na fábrica do Seixal da Siderurgia 
Nacional — Trata-se de uma cadeia tipo Dwight- 
“Lloyd com uma capacidade de aproximadamente 
250 000 tons/ano; 

b) Sinterização na fábrica de Vicominas em Vila 
Cova do Marão. Cuba fixa Greenawalt com uma 
capacidade diária de = 100 tons (considerando os 
concentrados ali produzidos); 

c) IPeletização usando dois discos de 45m de diã- 
metro em Vila Cova do Marão para a produção 
de peletes grangcold e com uma capacidade de 
60:70 tons/hora. 


Portanto uma capacidade de aglomeração insigni- 
ficante. 

No caso de Vila Cova do Marão a capacidade é 
excedente tendo em atenção à consumo do forno de 
redução. No entanto se a industrialização do País se 
vier a processar e o consumo de ferro fundido subir 
para valores próximos dos europeus a capacidade 
de aglomeração disponível poderá vir a ser aprovei- 
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No caso do Seixal a capacidade de aglomeração 
é insuficiente para alimentar o alto-forno. Daí que 
haja necessidade de carregar minérios em bocados 
que o País importa. 


5 — PERSPECTIVAS 


De momento estã por aproveitar o jazigo de Mon- 
corvo cujas reservas são vultuosas, Elas devem as- 
cender a 200 milhões de toneladas de concentrados 
com cerca de 64 % de Fe. Dado que se torna necessá- 
rio moer o minério bastante fino (100 % — 80 j) 
para o concentrar, o único processo de aglomeração 
aplicável será a peletização, De momento ainda se 
não conhece qual o tipo de peletes a fabricar: ácidas 
cu autofundentes, Tal só se tornará possível quando 
so [fizerem ensaios para se determinarem as caracte- 
rísticas de um e outro tipo e quando se fixar em defi- 
nitivo o lecal de implantação da instalação. Até lá no 
entanto parece-nos que há um longo caminho a per- 
correr, só não se compreendendo que a era das via- 
gens interplanetárias esse caminho não seja percor- 
rido em curto espaço de tempo. 


6 — INSTALAÇÕES LABORATORIAIS 


Do nosso conhecimento apenas existem na Fa- 
culdade de Engenharia da Universidade do Porto apa- 
relhos laboratoriais de aglomeração dimensionados 
para, a partir deles, se tirarem elementos suficientes 
para que da transposição para a escala industrial 
não surjam surpresas desagradáveis. 

A cuba de sinterização permite trabalhar com 
altura de carga até 50 cm e está equipada com dis- 
positivos de medição e registo de temperatura e de- 
pressão aos diferentes níveis de carga. Vários pro- 
cessos de mistura são disponíveis a fim de se estudar 
a sua influência na permeabilidade a frio e durante a 
aglomeração. A instalação é bastante flexível e per- 
mite obter depressões até 1350 mm de coluna de água. 

Os sínters são analisados visualmente e metalo- 
graficamente estando em curso de montagem uma 
instalação para determinação da reductibilidade, in- 
chaço e resistividade a frio e quente, 

Nesta instalação foram feitos ensaios de sinte- 
rização de minérios de ferro de Vila Cova com o fim 
de procurar melhorar a capacidade de instalação 
industrial. Mais recentemente têm sido realizados en- 
saios de aglomeração dos finos de minério de man- 
gamés e das escórias da metalurgia do ferro manga- 
més, com a finalidade de determinar a sua aptidão 
a aglomeração e a composição da carga em combus- 
tível, teor em humidade, etc. 

Sobre estes «sinters» realizam-se os ensalos ha- 
bituais já disponíveis pensando-se a curto prazo de- 
senvolver os emsaios de redutibilidade, inchaco e 
resistividade. 

Quanto à peletização dispõe-se de dois discos um 
de 20 cm de diâmetro e outro de 100 em de diâmetro, 
Este dispõe de 3 margens a fim de se lhe aumentar 
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a flexibilidade. Com efeito, quando foi construído 
admitia-se que o processo de flutuação viesse a im- 
por-se na concentração das hematites e que houvesse 
a necessidade de fabricar peletes envolvidas em finos 
de coque, Nesta instalação de dimensão piloto podem 
fabricar-se peletes de diversos calibres uma vez que é 
permissível alterar o débito da alimentação, o débito 
d: água, a velocidade e a inclinação do disco. 

As peletes são ensaiadas em verde à compressão, 
à queda e aos choques repetidos ao mesmo tempo que 
sº determina o teor em humidade. A sua cozedura 
se necessária, faz-se na cuba de sinterização por meio 
da gãs. As peletes cozidas são ensaiadas à compressão 
e observadas macro e microscopicamente. 

Nesta instalação foi recentemente realizado o 
ensaio de peletização do minério de Vila Cova pelo 
processo grangcold, com plemo êxito, A transposição 
do processo para a escala industrial não ofereceu 
quaisquer dificuldades nem houve grandes surpresas. 

Nesta mesma instalação se tem investigado o 
processo de fabrico das peletes pré-reduzidas por 
incorporação de finos de coque na mistura. Os resul- 
tados são de molde a que se prossiga, No entanto, 
torna-se necessário instalar um forno rotativo de 
cozedura e pré-redução, pois até aqui estas operações 


têm sido realizadas em cadinho sem possibilidade de 
controlo rigoroso de temperatura e de atmosfera ga- 
sosa O que provoca a formação de cachos, No entanto 
o grau de redução é muitíssimo elevado. 


7 — CONCLUSÕES 


Pretendemos apenas apresentar a aglomeração 
de uma maneira genérica uma vez que nos pareceu 
ser um congresso deste gênero pouco propício, para 
já, à aceitação de trabalhos de investigação numa 
faixa de conhecimentos necessariamente estreita. 

Foi também nosso propósito informar o que já 
existe, o que deveria existir e quais as perspectivas. 

Pelo que fica dito a peletização será o processo 
d: aglomeração que de acordo com o que se conhece 
virá a ser implantado no mundo, cada vez em maior 
escala. Será também o processo susceptivel de vir a 
ser instalado em Portugal para aproveitamento dos 
jazigos de ferro actualmente conhecidos. 

A aglomeração por sinterização terá maior possi- 
bilidade de se desenvolver entre nós, aplicada a finos 
de minérios de manganés e às escórias da metalurgia 
do ferro-manganés. 


NOTICIÁRIO 


e Ponte sobre o rio Guadiana em Vila Real de Santo António 


O projecto da ponte internacional sobre o rio Guadiana entre Vila Real de Santo António e Aymonte 
encontra-se concluído, tendo sido já entregue na Junta Autónoma de Estradas. 

A aprovação do projecto pelos Governos dos dois países terá lugar, conforme disposição do convénio 
luso-espanhol, em face da sua apreciação pela Comissão Técnica Mista Luso-Espanhola, cuja reunião já foi 


solicitada pela representação portuguesa, 


A presente obra é constituída, para além da ponte propriamente dita e que é encargo dos dois países, 


por um extenso viaduto na parte portuguesa de cerca de 500 metros, ao qual se seguem cerca de 6 quiló- 
metros de estrada com características de via rápida e que fará a ligação à E. N. 125 que actualmente costeia 
o Algarve. Simultaneamente serão edificadas as respectivas instalações alfandegárias. Os projectos das obras 
complementares encontram-se prontos e prestes a serem aprovados superiormente. 

Espera-se que esta obra seja iniciada no presente ano, prevendo-se a sua conclusão em 1980, satisfazendo 
assim uma das grandes aspirações da população algarvia, 

O projecto que vai agora ser entregue à Direccion General de Carreteras, da autoria do Prof. Edgar 
Cardoso, compreende o estudo de duas soluções estru turais, com vãos intermédios de 130 metros e extremos 
de 85 metros. 

A solução A caracteriza-se por a superestrutura ser constituída por 5 quadros ou pórticos fechados 
com 90 metros de envergadura total, em forma de trapézio, com a base maior para cima e prolongada em 
consola para um e outro lado, a fim de receber 6 tramos simplesmente apoiados, com 40 metros de vão. 

A solução B é constituída por uma superestrutura tipo viga continua, ou mais exactamente, pórtico 
múltiplo apoiado em livre dilatação condicionada no encontro do lado de Espanha e um pilar de transição 
do lado de Portugal e encastrado em 5 pilares gêmeos flexiveis. 

Para qualquer das soluções são apresentadas três hipóteses de fundação indirecta: poços cavados por 
havage, estacas verticais de betão armado, moldado no local, de grande diâmetro (2,00 m), estacas cravadas 
verticais inclinadas, com entubamento metálico depois cheio de betão armado, de médio diâmetro (1,00m). 

São os seguintes os custos orçamentados para as duas soluções estruturais: 

Solução A (fundações indirectas em poços), 298 000 000800 — alternativa AVI; fundações indirectas em 
estacas — grande diâmetro, 310 000 000$00 — alternativa AV2: (fundações indirectas em estacas — médio dia- 
metro, 316 000 000800. 

Solução B (fundações indirectas em poços), 296 000 000800 — alternativa BV1; (fundações indirectas em 
estacas — grande diâmetro), 312 000 000$00 — alternativa BV2: (fundações indirectas em estacas — médio diã- 
metro), 320 000 000800. 


(Informação de «O Primeiro de Janeiro», de 24-1-78) 
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Determinação de coeficientes de repartição sólido/líquido 
em ligas ternárias Fe-Mn-C 


RESUMO 


Usando a técnica de análise por microsonda, de- 
terminaram-se coeficientes de repartição sólido, líqui- 
do para algumas ligas da região rica em Fe do siste- 
ma Fe-Mn-C. O método utilizado envolveu a micro- 
análise de amostras temperadas a partir de tempera- 
turas de equilíbrio na região L+y. 


1 — INTRODUÇÃO 


O conhecimento dos valores experimentais dos 
coeficientes de repartição é muito importante em 
quaisquer estudos de solidificação. No caso dos sis- 
temas binários, esses valores podem ser obtidos a 
partir dos correspondentes diagramas de equilíbrio. 
Já para os sistemas ternários, os coeficientes de re- 
partição só podem ser calculados a partir dos dia- 
gramas de equilibrio quando estes incluem as linhas 
de conjugação («lie-lines:), o que geralmente não 
acontece visto que as técnicas clássicas empregadas 
no estudo de diagramas não permitem a obtenção 
dessas linhas [1], [2]. 

Uma das aplicações possíveis da análise por 
microsonda electrónica em estudos de diagramas de 
equilíbrio é precisamente a determinação das linhas 
de conjugação. Contudo existe muito pouca informa- 
ção publicada sobre o assunto, nomeadamente no que 
respeita ao equilíbrio sólido/liquido, onde o emprego 
da microsonda é presentemente a única técnica expe- 
rimental disponível. 

O estudo do equilíbrio sólido/liquido em ligas 
ternárias de Fe-X-C tem sido objecto de particular 
interesse, tendo sido já publicados dados relativos 
a coeficientes de repartição nos sistemas Fe-Cr-C 
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ABSTRACT 


Solid /liquid partition coefficients of some Fe-Mn- 
-€ allows, at the Fe — rich corner, have been deter- 
mined by means of electron probe microanalysis. The 
method used involved the microanalysis of specimens 
quenched from an equilibrium temperature in the 
L+y field. 


[3a 7] e Fe-Ni-C [8]. O presente trabalho destinou-se 
a iniciar a obtenção de dados similares para o siste- 
ma Fe-Mn-C. 

A técnica utilizada consistiu em determinar, atra- 
vés da microanálise, a composição das fases sólida 
e liquida, em equilibrio, presentes a uma temperatura 
situada na região (sólido + líquido), em amostras 
temperadas a partir dessa temperatura. A região 
estudada foi limitada à (líquido +-y) e os coeficientes 
de repartição do manganês e do carbono — Ko(Mn) 
e Ko(C), respectivamente — calculados a partir das 
composições assim obtidas. 


2 — EXPERIMENTAÇÃO 


21 — Ligas 


De uma série de ligas inicialmente preparadas, 
com composições destinadas a cobrir de uma maneira 
sistemática a região rica em Fe do diagrama Fe-Mn- 
-C, seleccionaram-se as de composição nominal indi- 
cadas na tabela I. 


(1) Assistente especialista do Laboratório de Fisica e Engenharia Nucleares, Sacavém. 
Comunicação apresentada ao Congresso da Ordem dos Engenheiros, Lisboa, Novembro 1977. 
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TABELA | — Composição das ligas 


N.º de Composição nominal 
referência (% em peso) 
1 0,35 € — 2 Mn 
2 1,090 — 2 Mn 
5) 1,600 — 2 Mn 


* Sic 0,02% 


As ligas foram preparadas em forno de indução 
d: 2Kg, em vácuo, a partir de ferro e manganês 
electrolíticos e de varões de grafite de alta pureza, 
e a partir delas obtidos varões com 16 mm de diã- 
metro dos quais se maquinaram os provetes necessá- 
rios. À composição final do material, na forma de 
provetes magquinados, encontra-se também na ta- 
bela T. 

A forma e dimensões dos provetes foram espe- 
cificamente condicionadas pela geometria do sistema 
montado para os ensaios com têmpera em bloco de 
cobre arrefecido, Por uma questão de uniformidade, 
foi utilizado também o mesmo tipo de provetes nos 
cnsaios com têmpera em água, 


2 2— Aparelhagem 


Utilizou-se uma montagem experimental especial- 
mente preparada para trabalhar em estudos de soli- 
dificação [7] e que consistia essencialmente num for- 
no de tubo vertical com possibilidade de circulação 
de argon e permitindo cfectuar têmperas directas. 

O forno, com tubo de mulite e clementos de aque- 
cimento de carboneto de silício, podia atingir tempe- 
raturas até 1600ºC, A temperatura avaliava-se por 
meio de um termopar associado a um controlador 
que permitia manter a zona quente a menos de 0,5º € 
do valor desejado, 

Quanto às têmperas a efectuar, foram utilizadas 
as duas técnicas designadas por A e E, a seguir 
descritas: 


TUBO DE ALUMINA 


TUBO PROTECTOR DO TERMOPAR E BUJÃO 


PROVETE 


FONTA DO TERMOPAR 


|) | 
27 | - 
E A ER 
fd ú mas 


ia ORIFÍCIO DE DESCA ISA 
Fig. 1 
Montagem do provete para têmpera em bloco de cobre 
arrefecido 


306 


Composição real (% em peso) * 


o Mn 


0,34 + 0,01 21 + 0,02 
0,96 + 0,01 2,07 + 0,02 
2,05 + 0,02 


1,55 + 0,02 
A — Têmpera em bloco de cobre arrefecido 


A têmpera efeciua-se sobre um bloco de cobre 
composto Ge duas partes articuladas e provido de 
refrigeração a água, situado abaixo da boca do forno. 
A parte inferior do bloco é mantida estática e a su- 
perior cai sobre a inferior logo que as primeiras 
gotas de metal nela se depositam. 

O provete a ensaiar e um tubo fino de alumina 
que vai actuar como protector do termopar e também 
como bujão («stcpper») são introduzidos num tubo 
de alumina mais largo provido de um orifício de des- 
carga, tal como se mostra na figura 1. 

Este conjunto permite colocar o provete no cen- 
tro da zona quente do forno, através de uma mon- 
tagem em que se fixa ao topo superior do tubo largo 
de alumina uma peca de latão refrigerada a água, 
obtendo-se assim uma mentagem rigida que, quando 
apoiada numa plataforma de fibrocimento colocada 
a uma dada altura fixa em relação ao tubo do forno, 
irá selar o extremo superior deste, Essa peça de 
latão dispõe também de uma tubagem que permite 
a circulação de argon purificado dentro do tubo que 
contém o proevete. O extremo inferior do tubo do 
forno é fechado por meio de uma outra peça metálica 
refrigerada a água e provida de uma tampa desli- 
zante. 

A têmpera efectua-se removendo essa tampa e 
levantando o bujão. A massa fundida ao cair acciona 
um sistema que força a parte superior do bloco de 
cobre a cair sobre o metal entretanto depositado na 
parte inferior, contribuindo assim para que o metal 
se espalhe mais rápida e uniformemente, o que con- 
duz a uma têmpera mais eficiente. 

Este tipo de têmpera é utilizado para tempera- 
turas de ensaio para as quais o material fundido é 
suficientemente fluido para permitir uma descarga 
eficiente. Das fases sólida e líquida em equilibrio à 
temperatura ensaiada resultam, na amostra obtida 
por têmpera, aquilo a que designaremos por «sólido» 
o «liquido». A estrutura deste «líquido; é bastante 
lina e permite uma Tácil determinação da sua compo- 
gição média, ao contrário do que acontece no «liquido» 
obtido por têmpera em água, 
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FIOS DE SUSPENSÃO DE PI 


TUBO PROTECTOR DO TERMOPAR 


TAMPA 


 . CADINHO DE ALUMINA 


Montagem do provete para têmpera em água 


PONTA DO TERMOPAR 


CIMENTO DE ALUMINA 


B — Têmpera em água 


O provete é introduzido num cadinho de alumina 
ao qual se adapta uma tampa e, centrado nesta, um 
pequeno tubo de alumina destinado a proteger a ponta 
do termopar (fig. 2). 

O cadinho é suspenso, por meio de fios de platina, 
do extremo superior de um tubo de alumina, idêntico 
ao tubo largo citado em A, e este adaptado ao tubo 
do forno tal como se referiu anteriormente, 

A têmpera efectua-se deixando-se cair o cadinho 
num recipiente contendo água, colocado por debaixo 
do forno. 

Este tipo de têmpera é normalmente associado 
a temperaturas de ensaio relativamente próximas do 
solidus, para as quais o material fundido já não 
apresenta grande fluidez, A estrutura do «liquido» 
obtido é geralmente grosseira, dificultando a deter- 
minação da sua correcta composição, mas por outro 
lado, uma vez que a amostra obtida por têmpera 
representa o total da massa fundida, é possível efec- 
tuar sobre ela a determinação das percentagens de 
fase «sólida» e «líquida» presentes, que serve de veri- 
ficação às composições de equilíbrio determinadas 
pela microanálise, 


2.3 — Ensaios 


Depois de efectuada a montagem correspondente 
ac tipo de têmpera a realizar, fazia-se passar uma 
corrente de argon no forno, durante cerca de 20 mi- 
nutos, antes de o ligar. O provete era então aquecido 
e fundido, mantendo-se neste estado durante 20 mi- 
nutos, ao fim dos quais a temperatura era descida 
lentamente até uma temperatura dentro da região 
(sólido + líquido), mantida nesse valor durante o 
periodo de tempo necessário para se atingir o equi- 
líbrio e depois se fazia a têmpera. O fluxo de argon 
era mantido durante todo o ensaio e reforçado um 
pouco, imediatamente antes da têmpera. 

Para cada uma das diferentes ligas estudadas 
houve necessidade de fazer sempre um primeiro en- 
saio a fim de determinar a temperatura «<liquidus> 
correspondente à respectiva composição, a partir da 
curva de arrefecimento. 

As temperaturas foram determinadas usando ter- 
mopares de Pt-Pt/13 5% Rh, calibrados a 900º € con- 
tra um termopar padrão, 

Para os tempos de 2 horas geralmente usados, 
as distâncias(*) correspondentes a equilíbrio contro- 
lado por um mecanismo de difusão são, para o Mn, 
da ordem dos 120 um na fase sólida e 12 mm na fase 
líquida, Os valores dos coeficientes de difusão do man- 
ganês usados neste cálculo — D=5x10-º cm?/seg na 
fase sólida y e D=5x10-" em2/seg na fase líquida — 
foram respectivamente calculado, para uma tempe- 
ratura e composição médias, a partir de uma ex- 
pressão dada por Wells e Mehl [9] e tomado como 
o valor típico para metais liquidos referido por 
Chalmers [10]. Previa-se então daqui que partículas 
de «sólido» de diâmetro igual ou inferior a 120 am 
e respectivo «líquido» adjacente estariam em equilí- 
brio em termos de 'Mn e por maioria de razão em C€, 
sendo portanto adequadas para a microanálise. 


2.4 — Amostras para a microanálise 


24.1 — Seleccionaram-se como padrões para as 
determinações de Mn e €, respectivamente: 


— Uma pequena porção de varão da liga 3 (com 
dimensões adequadas à montagem num porta- 
-padrões convencional) que foi encapsulada 
num tubo de sílica, em argon a aproximada- 
mente 1/3 de atmosfera, homogeneizada a 
1150" C durante duas semanas e temperada em 
água, 

— Uma amostra obtida de um ensaio em que se 
fez uma têmpera, em bloco de cobre arrefecido, 
de material dessa mesma liga 3, a partir de 
uma temperatura acima da de fusão. 


Porções de ambos os padrões foram novamente 
analisadas tendo-se confirmado que mantinham a 


1/2 
(*) X=2(Dt) onde X é a distância de penetração, num dado meio e no tempo t, de um elemento com coeficiente 


de difusão D nesse meio. 
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composição referida no quadro I e o resto montado, 
com pasta condutora de prata, em pequenos aros de 
aço de 5mm de diâmetro. Foram depois submetidos 
a um polimento convencional até pasta de diamante 
do 1/4 um, acabamento com alumina y, lavagem com 
água e álcool e secagem. 


O uso de padrões de ligas quer para o Mn quer 
para o €, em vez de padrões de elementos puros, 
justifica-se aqui pelo baixo teor destes elementos nas 
amostras a analisar. A utilização de padrões de com- 
posição próxima da do material a analisar é parti- 
cularmente importante no caso do carbono, dado que 
a grande diferença em número atômico entre este e 
os restantes elementos da matriz exigiria uma larga 
correcção, de eficiência duvidosa, se se usasse um 
padrão de elemento puro, 


24.2 —De cada ensaio com têmpera A resulta- 
ram várias porções de metal com espessuras da ordem 
dos décimos de milímetro; dessas porções selecciona- 
ram-se duas ou três que se montaram em conjunto 
a fim de constituirem uma amostra para a microaná- 
lise, De cada ensaio com têmpera E, a amostra para 
a microsonda obteve-se por montagem de metade da 
amostra resultante, cortada longitudinalmente, Como 
material de montagem utilizou-se normalmente ba- 
quelite condutora, de ferro. 


As amostras foram, em ambos os casos, polidas 
como se descreveu em relação aos padrões e contras- 
tadas depois ligeiramente com uma solução de me- 
tabissulfito de sódio a 35 %, a fim de revelar ag fases 
«sólida» e «líquida», Uma vez escolhidas as áreas 
e direcções adequadas à análise por microsonda, assi- 
nalaram-se com marcas de microdureza, e retirou-se 
o contraste pela repetição da parte final do polimento 
a partir da pasta de diamante de 1/4 qm. 


Nas figuras 3 e 4 apresentam-se exemplos de 
microestruturas típicas obtidas. 


Fig. 4 


Liga 3. Têmpera A, a partir de 1396.5º €, Contrastação 
com S,.0,Na, a 35 %. Ampliação 100 x 
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Fig. 4 


Liga 3. Têmpera BE, a partir de 1386º C. Contrastação com 
S.0,Na, a 35 %. Ampliação 100 x 


2.5 — Análise com microsonda electrônica 


Utilizou-se uma microsonda MICROSCAN Mk 5 
da Cambridge Instruments, com 2 espectrómetros dis- 
persores de comprimento de onda, trabalhando-se 
sempre com uma tensão de excitação de 20 KV e 
uma corrente de feixe que produzia uma corrente de 
0,5x10-s ampêres num padrão de ferro puro, 

Para ambos os elementos, Mn e €, a radiação 
medida foi a K. A cada elemento foi atribuído um 
dos dois espectrómetros, usando-se para o MnkKa um 
contador selado e um cristal analisador de fluoreto 
de lítio e para o CK um contador de fluxo, com uma 
mistura de argon/70 % metano, e um pseudo-cristal 
de estearato de bário, 

O pico Mn Ka foi sempre localizado na amostra 
a analisar e o CK num padrão de grafite. O fundo, 
quer para o manganês quer para o carbono, era me- 
dido num padrão de ferro, 

Para o CK usou-se o dispositivo discriminador 
a fim de excluir eventuais interferências por parte 
d: reflexões de ordem elevada dos elementos pesados, 
nomeadamente a de 2.º ordem do pico MnL4Ã, Usou-se 
também sempre o dedo de refrigeração («cold fin- 
ger») a fim de evitar, ou pelo menos reduzir, a con- 
taminação devida ao carbono, durante a análise. Ten- 
do-se verificado [11] que a taxa de contaminação 
quando o dedo de refrigeração está em operação é a 
mesma quer este esteja colocado directamente sobre 
a amostra quer não, foi esta última posição a usada 
durante todo o trabalho uma vez que permitia uma 
mais fácil utilização do sistema óptico da microsonda. 

A maior parte dos valores experimentais obtidos 
foi retirada de registos gráficos correspondentes a 
varrimentos pontuais («point scans» — varrimentos 
em que o feixe de electrões estã fixo e a amostra se 
meve mecanicamente ao longo de uma dada direc- 
cão), provenientes de um registador de 3 canetas 
que possibilitava o registo simultâneo das intensida- 
des Mn Ka e CK e ainda a monitoração da intensi- 
dade de corrente na amostra. 
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2.5.1 — Sobre os padrões referidos em 2.4 efec- 
tuaram-se varrimentos pontuais tendentes a avaliar 
da sua microheterogeneidade, tendo-se verificado ser 
esta desprezável., 

2.0.2 — Em cada amostra efectuaram-se varri- 
mentos pontuais através de três ou mais partículas 
da «sólido». Utilizou-se uma velocidade de 30 am/mia. 
e cada varrimento começava geralmente numa área 
da «líquido», atravessava diametralmente uma par- 
tícula de «sólido» e acabava em «liquido» novamente, 
sendo os comprimentos totais assim varridos da or- 
dem dos 500 um. As partículas analisadas tinham 
sempre dimensões da ordem dos 100-150 um, excepto 


quando se procurou deliberadamente verificar se a 
composição do «sólido» em manganês variava ou não 
com a dimensão das partículas, dentro de determi- 
nados limites. 

Cada varrimento nas amostras era seguido de 
outros menos extensos, à mesma velocidade, nos pa- 
drões de Fe, de Mn e de C. Nos casos em que a 
intensidade da corrente na amostra se mantinha par- 
ticularmente estável, este procedimento apenas se 
levava a cabo no início e no fim da série de varri- 
mentos efectuada numa mesma amostra. 

Nas figuras 5 a 7 apresentam-se alguns dos re- 
gistos típicos assim obtidos. 


Liga 3. Têmpera A, a partir de 1396.59 € 


Fig. 6 
Liga 3. Têmpera A, a partir de 1397º € 


Fig. 7 
Liga 3, Têmpera EB, a partir de 1346º O 
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A partir dos registos calcularam-se os valores de: 
(i) Coeficiente de repartição do manganês — Ko (Mn) 


Int. Mn Ka no «sólidos 
Ko (Mn) = 


Int. Mn Ka no «líquido» adjacente 


Entende-ze por intensidade do pico Mn Ka no 
«sólido», a intensidade constante medida no patamar 
correspondente ao atravessamento das partículas «só- 
lidas» deduzida do fundo medido no padrão de ferro 
e corrigida para o tempo morto. 


Entende-se por intensidade do pico Mn Ka no 
«líquido» a média dos valores da intensidade lidos 
na parte do registo correspondente ao «liquido» todos 
os 15 um, entre uma distância de 50 um e uma de 
cerca de 150 um para cada lado da partícula «sólida», 
média essa corrigida para o fundo e para o tempo 
morto. O facto de ze terem desprezado os valores 
incluídos nas referidas distâncias de 50 um mais pró- 
ximas do «sólido» justifica-se pela necessidade de 
evitar incluir na média calculada valores de inten- 
sidade que poderiam estar afectados por efeitos arti- 
ficiais de têmpera [5]. 


(ii) % Mn no «sólido» e no «liquido» 


Nas amostras provenientes de têmpera A onde 
a segregação no «liquido: é minima, a média das 
intensidades não corrigidas de Mn Ka no «líquidos 
calcula-se facilmente e dá valores concordantes qual- 
quer que seja a direcção em que o «sólido» é atra- 
vessado. Por outro lado, a intensidade não corrigida 
de Mn Ka nos patamares correspondentes a parti- 
culas «sólidas» mantêm-se constante, não variando 
com o tamanho destas até pelo menos 200-250 um. 
Consequentemente, esses valores de Mn Ka, depois de 
corrigidos para o fundo e para o tempo morto, foram 
usados no cálculo da percentagem de Mn no «liquido» 
e no «sólido», assumindo-se uma relação linear entre 
a intensidade da radiação X e a concentração do ele- 
mento e usando-se como padrão de Mn o referido 
em 2.4. 


Nas amostras provenientes de têmpera E, quer 
a pequena percentagem de «líquido: presente quer 
a estrutura dendritica grosseira nele desenvolvida con- 
tribuem para tornar difícil a obtenção de valores 
médios correctos de Mn Ka no «líquido», quando se 
utiliza um varrimento pontual. Nestes casos, a inten- 
sidade Mn Ka no «liquido» usada quer na determina- 
ção da % de Mn quer na determinação de Ko (Mn) 
foi obtida através de medições associadas a varrimen- 
tos de área («rasters:-varrimentos feitos com a amos- 
tra fixa e o feixe de electrões deflectido magmnetica- 
mente por forma a cobrir uma dada área). Cobriam-se 
então áreas da ordem dos 50 a 100 (um): associadas 
a contagens por um tempo fixo de 40 segundos e 
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efectuavam-se também varrimentos de áreas equiva- 
lentes no padrão de ferro, a fim de determinar o 
fundo correspondente. 


(ii) Coeficiente de repartição do carbono-Ko (C) 


o? 


Determinou-se de maneira análoga à referida 
para Ko(Mn), mas com as intensidades CK apenas 
corrigidas para o fundo. 


(iv) %C no «sólido» e no «liquido» 


Foi calculada a partir das intensidades CK usa- 
das para o cálculo de Ko(lC), assumindo uma relação 
linear entre a intensidade e a concentração e usando 
como padrão de carbono o referido em 2.4, 


3 — RESULTADOS 


Na tabela II sumarizam-se os ensaios efectuados 
e respectivos resultados. As composições que aí se 
referem para as ligas são as composições nominais, 
podendo a correspondência com as composições reais 
ser obtida a partir da tabela 1. 

Note-se a existência de um ensaio com têmpera 
em água (liga 2, temperatura 1442,5º C) diferente 
dos outros do mesmo tipo na medida em que corres- 
ponde a uma temperatura relativamente próxima do 
liquidus. A proporção de «líquido» presente na amos- 
tra resultante é portanto muito mais elevada do que 
nas restantes amostras que sofreram o mesmo tipo 
de têmpera e, consequentemente, as considerações 
relativas à utilização de varrimentos de área feitas 
em 2.0.2 não se aplicam a este caso, A realização 
deste ensaio foi devida ao insucesso de várias tenta- 
tivas anteriores, de têmpera em bloco de cobre arre- 
cido, para a mesma liga. 

Relativamente à liga 1, os dados obtidos não 
permitiram estabelecer conclusivamente se foi de 
facto ensaiada na região (líquido + y). Esta liga não 
foi portanto considerada nas figuras 8 e 10, embora 


os valares resultantes da sua microanálise constem 
da tabela TI. 


3.1 -—Os resultados da microanálise efectuada 
sobre as amestras estão condensados na tabela II sob 
a forma de concentrações em Mn e € do «sólido» e 
do «liquido» em equilíbrio a uma dada temperatura, 

3.2 — Os pontos correspondentes às composições 
eliquido» e «sólido» para as temperaturas ensaiadas 
estão representados na figura 8, juntamente com 
valores liquidus e solidus dos diagramas binários 
FesC [12] e Fe-Mn [13]. 
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P liquidus 


Fo | Do. RR? 


V solidus 


da A o Mn % 


Fig. & 
Aspecto das superfícies y — liquidus e y — solidus na região rica em Fe do sistema Fe-Mn-C. Composições 
em % em peso e temperaturas em ºC 


Não se representaram nesta figura os valores velocidade de arrefecimento demasiado grande (30º C/ 
liquidus resultantes dos ensaios prévios de análise “min.) quando comparadas com as condições normal- 
térmica. A montagem experimental utilizada nesses mente usadas numa análise térmica convencional. 
ersaios não era propriamente a mais indicada visto Da mesma maneira, não havia facilidades de agita- 
que as cargas eram demasiado pequenas (l2g) ea ção do banho fundido, nem se tomaram precauções 
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Linhas de conjugação da região L 


KLC) KCMn) 


4 ó 75 
cid — Q 75 
is— -— .T 
O) 
45 - 17 
«42 - O — 76 
Õ3-— E —— «79 
Eos O — 78 


Fe | 2 3 Mn % 
Fig. 9 


Ko (Mn) e Ko (C) representados em pontos correspondentes a 50% de sólido, 
Composições em % em peso 


di g “4 


Fig. 10 
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Mn % 


y do sistema Fe-Mn-C. Composições em % em peso e temperaturas em ºC 


especiais para evitar a volatilização do Mn ou se 
dispôs de uma análise quimica do banho tal qual, 
não havendo portanto dados quanto à composição 
homogénea exacta à qual a temperatura liquidus de- 
terminada correspondia. 

SS — Os coeficientes de repartição da tabela II 
encontram-se representados na fig. 9, juntamente com 
valores de Ko(C) e Ko(Mn) correspondentes aos sis- 
temas binários Fe-C e Fe-Mn. A mesma informação 
é apresentada, sob outra forma, na fig. 10. 


4 — DISCUSSÃO 


4.1 — Quando se comparam as características ge- 
rais dos registos gráficos obtidos verifica-se que, nas 
amostras provenientes de têmpera em água (fig. 7) 
a segregação no «líquido» é muito mais acentuada 
do que nas provenientes de têmpera em bloco de 
cobre arrefecido (figs. 5 e 6). 


Outra das diferenças manifestada pelas amos- 
tras provenientes dos dois diferentes tipos de têmpera 
reside no gradiente de concentração do soluto no «só- 
lido». Parte deste gradiente pode ser atribuído à fun- 
ção de sonda [14] e ele pode ser também influenciado 
pelo relevo próprio da amostra, pelas fases localizadas 
abaixo da superfície e, pelo menos no caso do car- 
bono, pela constante de tempo utilizada. Quando se 
tomam em consideração estes factores, pode consi- 
derar-se que a concentração em Mn no «sólido» não 
apresenta gradiente real e, consequentemente, as par- 
tículas de «sólido» estão em equilíbrio com o «líqui- 
do». O gradiente em carbono, quando existe, será 
então devido à ineficiência da têmpera e tal como 
previsto é mais acentuado em amostras provenientes 
de têmpera em água. 

Quanto à possivel existência de efeitos artificiais 
d> têmpera [5] o registo da fig. 6 não apresenta 
evidência deles, enquanto os picos mais elevados que 
delimitam o «sólido» na fig. 7 tanto podem ser atri- 
buídos a esses efeitos como não, O mesmo se poderá 
dizer dos dois picos mais acentuados que se situam 
de um dos lados do «sólido» no registo da fig. 5 e, 
genericamente, de todos os outros registos obtidos. 
Constata-se portanto que este tipo de efeitos não 
foi identificado de maneira conclusiva no presente 
trabalho, o que pode ter resultado não só de têmpe- 
ras mais eficientes mas também de uma menor sen- 
sibilidade na detecção desses efeitos quando se utiliza 
um varrimento pontual em vez de um varrimento de 
linha como em [5]. 


4.2 — Os poucos valores liquidus e solidus obtidos 
apenas permitem, com o auxilio dos correspondentes 
valores dos sistemas binários, o estabelecimento da 
tendência geral das linhas isotérmicas nas superfícies 
y-liquidus e y-solidus, próximo das composições estu- 
dadas (fig. 8). 

Para o sistema Fe-Mn+-C, os únicos valores liqui- 
dus publicados correspondem a concentrações relati- 
vamente altas de soluto (Mn > 11 %, em peso) [15]. 
Uma vez que não existe mais informação sobre o 
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assunto, não há possibilidade de se fazerem compara- 
ções com os valores agora obtidos, podendo apenas 
confirmar-se que a tendência geral das isotérmicas 
da fig. 8 está de acordo com o trabalho de Isobe [15]. 
Existe contudo um valor (1442,5º C) que não se ajusta 
a essa tendência, o que sugere que pode ter havido 
um erro na determinação da composição do «líquido» 
correspondente a esse ponto, 

Em relação à superfície y-solidus, não hã infor- 
mação publicada com a qual se possa estabelecer 
comparação. 

44 -—Da comparação entre os coeficientes de 
repartição correspondentes aos dois primeiros ensaios 
referidos na tabela TI, efectuados aproximadamente 
à mesma temperatura mas com diferentes tempos, 
pode concluir-se que o tempo de 2 horas foi já sufi- 
ciente para se estabelecer o equilibrio. Pode portanto 
inferir-se que em todos os outros ensaios este foi 
também atingido, pelo menos nas partículas «sólidas» 
com dimensões dentro dos limites considerados. 


Em geral não se encontravam numa mesma amos- 
tra partículas «sólidas» abarcando uma grande gama 
d: tamanhos mas, nos casos em que se efectuaram 
varrimentos sobre partículas de dimensões variando 
entre 75um e 250 yum, não se detectou nenhum de- 
créscimo sistemático na concentração em Mn com o 
acréscimo de dimensão. Uma vez que a difusão do 
carbono é mais rápida que a do manganês, a cons- 
tância do teor em Mn no «sólido» pode ser tomada 
também como prova de que o equilíbrio foi de facto 
atingido, Isto confirma o que se referiu anteriormente 
e mostra que as distâncias calculadas em 2.3 foram 
subestimadas, 

Da fig. 9 conclui-se que Ko(Mn) (para uma per- 
centagem fixa de Mn de cerca de 2 % em peso) não 
varia significativamente com o teor em carbono. Re- 
portando-nos às conclusões de Philibert de que o car- 
bono aumenta o índice de segregação do elemento X 
em ligas Fe-X-C [16], confirmadas pelo trabalho de 
Doherty e Melford em ligas Fe-Cr-C [3] a explicação 
ai sugerida segundo a qual esse efeito era provavel- 
mente devido à diminuição de Ko(X) não parece por- 
tanto ser válida no caso de X=Mn. 

Quanto a Ko(C) aumenta com a presença do 
manganês e, se se excluir o valor anómalo de 0,53 
(correspondente à temperatura de 1442,5” C anterior- 
mente referida), parece ter tendência a aumentar 
também ligeiramente com a percentagem de carbono. 
Contudo este aumento pode resultar de um efeito real 
de temperatura em Ko, o que não seria de estranhar 
dado que as superfícies y-liquidus e y-solidus não são 
planas. 


CONCLUSÕES 


Confirmou-se que a técnica da microanálise é de 
grande importância na determinação experimental 
dos coeficientes de repartição entre sólido e líquido. 
Utilizando essa técnica e dentro de uma gama muito 
limitada de composições na região (L + y) do sistema 
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Ensaios 


TABELA II 


Ko(C) 


0,45 
0,45 
0,46 
0,46 
0,46 
0,53 
0,42 
0,39 


Carbono (% em peso) 


2,03 
1,95 
2,26 
2,06 
2,31 
1,07 
1,93 
0,83 


«Liquidos 


0,92 
0,88 
1,04 
0,94 
1,06 
0,57 
0,64 
0,32 


«Sólidos 


KoiMn 


0,78 
0,77 
0,75 
0,79 
0,75 
0,79 
0,76 
0,78 


em peso) 


Manganês (% 


 ——— 


2,15 
212 
2.18 
212 
2,36 
2.02 
2.25 


«Liquido 


«Sólidos 
1,68 
1,64 
1,64 
1,68 
1,79 
1,61 
1,72 
1,86 
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Fe-Mn-C, verificou-se que, em relação aos sistemas 
binários Fe-C e Fe-Mn: 


— O coeficiente de repartição do manganês não 
é afectado significativamente pelo teor em car- 
bono; 

— (O) coeficiente de repartição do carbono aumen- 
ta com a presença do manganês, 
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Solução do trânsito de potência 
pelo método do desacoplamento 


J. M. FERREIRA DE JESUS (1) 
J. P. SUCENA PAIVA (1) 


RESUMO 


O cáleulo do trânsito (ou fluxo) de potência numa rede de energia eléctrica, é presentemente efectuado 
em computador digital por meio de diversos métodos numéricos, dentre os quais o mais sofisticado é o 
método de Newton-Raphson e suas variantes, 

No artigo os autores descrevem um programa de cálculo baseado numa destas variantes, o método 
do desacoplamento, inicialmente proposto por Stott [10]. 

É dada especial ênfase às técnicas de programação com esparsidade utilizadas. Apresentam-se os resul- 
tudos numéricos obtidos para diversas redes, comparando-as com os obtidos pelo método de Newton-Raphson 
clássico. 


SUMMARY 


The power flow solution in an electric energy system is currently executed by digital computer using 
several numerical methods, of which the most powerful is the Newton-Raphson and its variations. 
In this paper the authors describe a computer programme based upon one of these variations, the 
decoupling method, initialhy suggested by Stott [10]. 
The sparsity programming thecniques are specially emphazised. Numerical results for several real 
systems are presented, and compared with those obtained by the classic Newton — Raphson method, 


1 — INTRODUÇÃO 


O cálculo de uma rede de energia eléctrica em regime estacionário, geralmente designado por trânsito 
(ou fluxo) de potência, é presentemente efectuado em computador digital. Antes do advento deste, e ainda 
durante algum tempo depois, a solução do trânsito de potência era obtida por meio do analisador de redes. 

As primeiras soluções por computador digital apareceram na literatura em 1956 [1] e 1957 [2] e 
empregavam o método dos deslocamentos sucessivos, geralmente conhecido por Gauss-Seidel, Este método 
iterativo, que utilizava a representação da rede pela matriz das admitâncias nodais, em breve se tornou 
inadequado devido à sua reduzida velocidade de convergência. 

O método de Newton-Raphson, inicialmente proposto por Van Ness [8], apresenta características de 
convergência largamente superiores. No início da sua aplicação à solução do trânsito de potência revelou-se 
pouco competitivo devido às suas exigências no que concerne à memória ocupada, Com o aperfeiçoamento 
por Tinney e outros [4] — [7] dos métodos de programação com esparsidade, o tempo de cálculo do método 
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de Newton-Raphson foi substancialmente reduzido e as necessidades de memória diminuidas, tornando este 
método o mais sofisticado que o engenheiro de sistemas tem presentemente à sua disposição para a solução 
do fluxo de potência. 

Devido ao aumento continuo da dimensão das redes de energia a calcular e ao uso cada vez mais 
generalizado dos programas de trânsito de potência, a investigação prosseguiu com a vista a tornar o método 
de Newton-Raphson ainda mais eficiente. Uma característica essencial de qualquer rede de energia é a forte 
interdependência entre o trânsito de potências activas e os argumentos das tensões, por um lado, e entre o 
trânsito de potências reactivas e os módulos da tensão, por outro. Em contrapartida o acoplamento entre potên- 
cia activa e módulo da tensão, e potência reactiva e argumento da tensão é relativamente pouco importante. 
Recentemente Stott [9] [10] tirando partido destas propriedades propôs o método do desacoplamento no qual 
os trânsitos de potência activa e reactiva são resolvidos separadamente, 

Em relação ao método clássico de Newton-Raphson, o método do desacoplamento, que se pode consi- 
derar uma variante daquele, necessita de um número em geral mais elevado de iterações até atingir a con- 
vergência. Contudo o tempo de cálculo por iteração é substancialmente (cerca de 5 vezes) menor, o que 
resulta numa apreciável redução do tempo de cálculo. As exigências de memória são também consideravel- 
mente (cerca de 40%) inferiores às do método de Newton-Raphson. 

Nesta comunicação descreve-se um programa de cálculo executado no Centro de Electrotecnia das Uni- 
versidades de Lisboa, que utiliza o método do desacoplamento proposto por Stott, com algumas variantes. 

Dado que é difícil encontrarem-se na literatura detalhes das técnicas de programação com esparsi- 
dade, que são a chave do êxito dos métodos de cálculo expostos, pareceu aos autores que teria interesse desen- 
volver com algum pormenor este aspecto, | 

A utilidade dos métodos de programação com esparsidade não se restringe aliás à solução do trânsito 
de potência, encontrando aplicação nos mais diversos ramos da engenharia. 

O programa foi aplicado ao cálculo diversas redes, tendo-se sempre obtido a convergência com rapidez 
e segurança. No corpo do artigo dá-se conta dos resultados obtidos, que são comparados com os do método 
de Newton-Raphson clássico. 


2 — METODO DO DESACOPLAMENTO 
As equações do fluxo de potência numa rede de energia em regime estacionário são [15] 


n 


P,-10,-D 9, V, V,=0 (1) 


Representando as tensões sob forma polar, as n equações complexas [1] transformam-se em 2n equa- 
ções reais 


F=U, » U, G,, cos (5, — 85) + B, sin (5, — 84) = (2.8) 
e. | 
| 

o — 1, > U, G, sin (4, -—8,) —- B, cos (8, — 8,) = 0 (2.b) 
j=i | 


| = ; poe À 
A cada barramento da rede corresponde um par de equações [2]. Há porém que considerar que numa 
rede de energia se definem três tipos de barramentos 


— Barramentos tipo PQ, para os quais são conhecidas as potências injectadas P e Q, e desconheci- 
dos o módulo e argumento da tensão U es.. 

— Barramentos tipo PU, para os quais são conhecidos a potência activa injectada P. e o módulo da 
tensão U, e desconhecidos argumento da tensão ô, e a potência reactiva Q.. Estes nós podem ser 
transformados em nós tipo PQ se durante o cálculo se verificar que Q, excede os limites do gera- 
dor ligado ao nó i, Qin € o 
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— Barramentos de balanço ou flutuantes, para os quais são conhecidos o módulo e argumento da ten- 
são e desconhecidas as potências activa e reactiva, Em geral numa rede existe apenas um nó deste 
tipo, podendo no entanto, haver mais do que um, 


A aplicação do método de Newton-Raphson à solução das egs. [2] conduz a [15] 


SP H N às 


(3) 


| 


AQ - E AU/U 


Os erros de fecho AP, e 4Q, são respectivamente os valores dos primeiros membros das equações (2a) 
e (2b). 


Os elementos do Jacobiano são dadas pelas expressões seguintes [15]: a 
Para 1) 
H,,=L,=U, U,[ G,sin (8-8,))-B, cos (8/-3,)] (4a) 
Ny; =J,;=U, U, [ G, cos (3,-8,)) +B, sin (8, —8,)] (4b) 
Para i=j 
H,,=-Q,-B,, U, (5a) 
Li= Q4-ByU (5b) 
Ni= P+G, U; (De) 
= R-GU (5d) 


Se estas equações forem satisfeitas com determinado conjunto de valores dos argumentos 5, e módu- 
los U, das tensões nodais, então AP = 4Q, =0 e o problema estã resolvido. Caso contrário a equação matri- 
cial [8] permite calcular correcções 45, e AU, a adicionar aos valores de 5, e U,. Desta operação resultam 
melhores aproximações dos valores correctos. O processo iterativo dá-se por terminado quando AP e AQ 
forem inferiores a uma dada tolerância especificada (por exemplo 0,01 MW/MVAR). 

Os erros de fecho para os barramentos de balanço não são incluídos na equação [8] nem AQ, e àU, 
para os nós do tipo PU. A ordem do sistema de equações o resolver é por conseguinte (Dt 2 Dpq)- Se os 
nós i e j não estão directamente ligados então o elemento Y;, da matriz das admitâncias nodais é nulo,, 
o mesmo sucedendo aos elementos correspondentes das submatrizes [HJ], [NJ], [J] e [L]. O Jacobiano é por 
conseguinte uma matriz muito esparsa com simetria topológica mas não numérica. 

O método de Newton-Raphson, na versão polar apresentada converge quase sempre em 2 a 5 itera- 
ções, independentemente do número de barramentos da rede, a partir de uma partida V,=1,0 L Eid PD. U. 
Se houver mudança de tipo de nós (PU para PQ ou vice-versa), devido à necessidade de satisfazer limites 
rígidos de potência reactiva ou de tensão, ou se proceder à actuação automática da relação de transfor- 
mação dos transformadores para regular a tensão em determinados barramentos o número de iterações 
aumenta consideravelmente. De salientar que o Jacobiano tem de ser calculado e triangularizado em cada 
iteração, o que representa uma percentagem substancial do tempo total de cálculo por iteração. 

O método de desacoplamento é uma variante do método de Newton-Raphson no qual se desprezam as 
submatrizes [N] e [J] do Jacobiano, que são a representação matemática do fraco acoplamento [AP] e 
[4U/U], e [AQ] e [45], respectivamente. Obtêm-se assim equações separadas para os fluxos de potência 
activa e reactiva 


[4P] = [H] [ 48 ] (6) 
[4Q] = [L] [AU/U] (7) 


Esta simplificação, embora permita uma redução substancial da memória exigida ao computador, 
resulta uma degradação das características de convergência do método de Newton-Raphson, pelo que o 
número de iterações até se atingir a convergência aumenta. A redução total do tempo de cálculo é pouco 
significativa. 
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O método só se torna verdadeiramente eficiente se a submatrizes [H] e [L] forem constantes, o que 
constitui a característica mais saliente do método proposto por Stott [10]. No Anexo 7.2 apresenta-se a 
formulação deste método, que leva até às últimas consequências o princípio do desacoplamento. As equa- 
ções [6] e [7] transformam-se em 


[AP/U] = [B'] [48] (8) 
[49/U] = [B"] [AU] (9) 
onde 
, 1 
R==—— HSM (10) 
ij 
n 
MR > a 11 
5 x, (11) 
j=1 1 
Ji 
bj = — B (12) 


As matrizes [B] e [B”] representam aproximações constantes aos declives dos hiperplanos tangentes 
às funções AP/U e AQ/T, respectivamente, São matrizes reais, esparsas e simétricas, com elementos cons 
tantes. 

Para a solução das equações (8) e (9), após uma comparação dos possíveis métodos, foi decidido 
utilizar a bi-factorização, inicialmente proposta por Zollenkopf [17]. A fim de minimizar a criação de novos 
elementos não nulos que a factorização acarreta, procede-se previamente a uma ordenação dinâmica das 
colunas de [B'] e [B']. 

Dado que [B'] e [B'"] tem elementos constantes a bi-factorização é efectuada apenas uma vez antes 
do início do processo iterativo. 

As equações (8) e (9) podem ser resolvidas simultaneamente em cada iteração, sendo contudo muito 
mais eficiente resolver primeiro a equação (8) para obter [45] e corrigir os valores de [8] antes de cal- 
cular [AQ/U] e resolver a equação (9) para obter [AU], Este método conduz a valores muito aproximados 
da solução apenas após uma ou duas iterações. As iterações subsquentes permitem refinar estes valores, 
sendo no entanto perdida a convergência quadrática do método de Newton-Raphson. 

Em comparação com o método de Newton-Raphson, cada iteração do método do desacoplamento é 
muito mais rápida embora o número total de iterações até se atingir a convergência seja, para a mesma 
tolerância dos erros de fecho, em geral mais elevado. O tempo total de cálculo vem consideravelmente redu- 
zido, bem como a memória ocupada. 


3 — DESCRIÇÃO DO PROGRAMA DE CÁLCULO 


O programa desenvolvido é composto por sete subrotinas de cálculo e três subrotinas entrada-saída. 
Descrevem-se em seguida as subrotinas de cálculo, dando especial atenção às técnicas de programação com 
esparsidade utilizadas. 


3.1 — Armazenamento da matriz de admitâncias nodais 


Sendo a matriz de admitâncias esparsa e simétrica, são apenas armazenados os termos não nulos 
correspondentes ao triângulo superior. Os elementos diagonais da matriz de admitâncias são armazenados 
num vector DIAG, 

Os elementos não diagonais do triângulo superior são armazenados num vector VALOR. Para detectar 
a posição destes elementos são criados dois vectores LINHA e ICOL, 

Tem-se assim duas tabelas: na primeira armazenam-se os elementos diagonais da matriz de admitân- 
cias. O índice do vector DIAG é o número do nó de que se pretende o elemento diagonal; na segunda tabela 
o índice de VALOR, LINHA e ICOL é um índice mudo. Em LINHA é armazenado o índice da linha a 
que o elemento VALOR pertence, sendo armazenado em ICOL o índice da coluna correspondente. 
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A im A 


Exemplificando, suponha-se que se pretende armazenar a matriz: 


[A] = 2 —1 —2 0 
=ri 3 O — 4 

— 2 O) 5 — 5 

0 —d —5 6 


As tabelas ficarão preenchidas como a seguir se indica; 


TABELA 1 
I | 1. o MS á 
DIAG (1) | 2 3 5 6 
TABELA 2 
K | 1 2 3 4 
VALOR (K) | st. cmi — 4 —5 
LNHA (K) | 1 1 
ICOL (K) | 2 3 4 4 


Os elementos da matriz das admitâncias têm interesse para a determinação de [5P/U] e [4Q/U], que 
ocupa uma parcela importante do tempo de cálculo por iteração. A tabela 2 é percorrida sequencialmente, 
desde K=1 até ao valor máximo, sendo detectados os elementos G,, e B, e calculados os termos U, U, 
G,eU,U, B,,. Estes termos podem ser usados directamente para o cálculo AP,, aP,, sq, e 4Q, (equações 
(2a) e (2b) por soma com os valores acumulados. 

As funções trignométricas sin (8, — 8,) = — sin (58, — 8,) e cos (5, — 8,) = cos (ô, — 3,) podem, para 
para tornar o cálculo mais rápido, ser aproximados para dois e três termos da expansão em série, respec- 
tivamente, sem perda significativa de precisão. 


3.2 — Armazenamento das matrizes [B'] e [B'] 


As matrizes [B'] e [B'] são simétricas e têm esparsidade idêntica à matriz de admitâncias, Dado 
que a factorização acarreta em geral a criação de novos elementos, [B'] e [B"] devem ser armazenados 
de forma a permitir o processamento dos novos elementos com eficiência. 


O método de armazenamento de [B'] e [B'] baseia-se nas duas tabelas a seguir indicadas: 


TABELA 1 


DIAG1i — valor numérico do elemento diagonal de [B'] 

DIAG2 — valor numérico do elemento diagonal de [B”] 

ILOC —localização do primeiro elemento não diagonal da coluna correspondente na tabela 2 
NANU número de elementos não nulos não diagonais na coluna correspondente 


TABELA 2 


VALORi — Valor numérico do elemento não diagonal de [B'] 

VALOR2 — Valor numérico do elemento não diagonal de [B"] 

ILINHA — indice da linha do elemento não diagonal de [B'] e [B'"] 

IPROX Localização do próximo elemento não nulo na coluna; se não existir próximo elemento 
é colocado um zero na posição correspondente ao último elemento da coluna, 
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A matriz [A] usada previamente será, de acordo com este esquema armazenado do modo seguinte: 


TABELA 1 
| Í 2 3 4 
DIAG (1) | 2 E 5 6 
ILOC (1) | 1 3 5 7 
NANU (1) | 92 2 2 2 
TABELA 2 
K | 1 2 3 4 5 6 7 8 9 10 11 
VALOR (K) | anel a A —4 saiá E: = 4 —5 0 
ILINHA (K) | 2 3 1 4 1 4 2 3 0 0 
IPROX (K) | 2 0 4 0 6 0 8 0 10 11 


Para poupar memória, os vectores auxiliares ILOC, NANU, ILINHA e IPROX são comuns às duas 
matrizes [B] e [B"], o que levanta algumas dificuldades no tratamento dos nós tipo PU. Com efeito na 
matriz [B"] não existem as linhas e colunas correspondentes aos nós deste tipo. Na ordenação e factori- 
zação os elementos destas linhas e colunas não são por conseguinte processados. Contudo se houver infracção 
dos limites de potência reactiva, o nó sofre um tratamento especial (ver secção 3.6), para o qual é neces- 
sário conhecer os valores destes elementos, 

Nestas condições é necessário criar mais dois vectores LOC e NPROX pertencentes respectivamente 
às tabelas 1 e 2. LOC indica a posição do primeiro elemento não diagonal da coluna correspondente ao nó 
PU na tabela 2; NPROX indica a posição do próximo elemento não diagonal da mesma coluna. 

A primeira vista a memória necessária para o armazenamento de [B'] e [B"] poderá parecer excessiva, 
nomeadamente o armazenamento de todos os elementos e não apenas os triângulos superiores, como no caso 
da matriz de admitâncias. Tal é contudo aparente, visto a informação contida nas duas tabelas criadas 
ser toda necessária a reconstrução de [B'] e [B”] e para a realização da factorização em qualquer sequên- 
cia de eliminação. Além disso durante a factorização são geralmente criados novos elementos não nulos, 
que são armazenados nas posições livres criadas durante a ordenação. Estas correspondem aos elementos do 
triângulo superior que no decorrer do processo forem eliminados. 


3.3 — Bifactorização 


No programa que ora se descreve utilizou-se a bifactorização das matrizes [B'] e [B'] na solução 
das equações (8) e (9). 

O método, desenvolvido por Zollenkopf [17], consiste na determinação de 2n matrizes factores (para 
uma matriz n x n) de tal modo que o produto destas 2n matrizes satisfaça a relação: 


[1,72] [LED + [1.27 [DI [A] [RP] [6º]: [Rº-D] [RP] = [5] 
(13) 
onde [A] —matriz original 
[L] — matrizes factores esquerdas 
[R] — matrizes factores direitas 


[1] — matriz identidade 


Pré-multiplicando sequencialmente a equação (13) pelos inversos de [4] , id à dis [1,27] e 


[4] + Vem 
[4] [RD] [Rº] [R (0-5) [Rº]= [9] [97 [=D [e (14) 
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Pós-multiplicando sequencialmente (14) por [L”]: [pg se ED] [4,8] VOS 


Ea Te PR TR e E] EC (15) 


-— 
Finalmente pré-multiplicando (15) por [A] vem: 


[R2] [8] E [RD] [R?] [1.8] [ES-D7 + [1.2] [1.27] = [A] (16) 


As matrizes factores permitem asim determinar a inversa de [A] a partir da expressão (16). 
A solução do sistema de equações 
[A] [X] = [B] (17) 


é por conseguinte dada por 


Pe Le a RI RS] E] Te] E] Da] (18) 


A dedução das expressões das matrizes factores não se apresenta por falta de espaço. Verifica-se [17] 
que 


Eu] mi 1 ),, 
] : 
13 (19) 
L, 1 
na da 1 
onde 
(k—1) 
E x 
jm (K=—1) na 
a. 
[5 
(20) 
| — E 
ij (K—1) 
15 
e 
KH) = 
[Rr | 1 
Fa 1 F F,n 
1 | (21) 
1 
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onde 


(K=-l) 
a 
É ESSE E (22) 
Já (K=1) ai 
ji 
O gas la fó la de s 
Os elementos d que figuram nas equações (20) e (22) podem ser obtidos pela fórmula de re 
corrência 
(K--1) (K —1) 
(E) E) CU EL ij 
a, =8a,) — 0 (23) 
à (K—1) 15 
8] 
(K-l1) (K—l) 
be são os elementos originais da matriz [A]. Se esta for simétrica, tem-se a,, =a,, , sendo 
por consequência r,, ep e, Pode concluir-se portanto que, àparte o elemento diagonal, [R] é a trans- 


posta de [Lj]. 


Da equação (21) pode verificar-se que a matriz [Rº | é a matriz identidade. A equação (18) simpli- 
fica-se portanto 


[x] = [RR] Escd Rs [nr] [1.872] Eds ES [2] [ut] [g| (24) 


A solução de um sistema de equações lineares por meio da bifactorização é basicamente equivalente 
a uma eliminação de Gauss-Jordan. Nas matrizes factores [L] e [R] está contida toda a informação rela- 
tiva às operações necessárias para efectuar esta eliminação. Do ponto de vista de cálculo a bifactorização 
é porém muito mais eficiente pois as matrizes [L] e [R] são calculadas uma só vez antes do início do 
processo iterativo e armazenadas apenas os elementos não nulos (com a excepção dos elementos diagonais 
unitários). Soluções repetitivas do sistema de equações (17) com vectores de termos independentes [B] dife- 
rentes são levados a cabo com extrema rapidez por utilização da equação (24). O número de posições de 
memória necessárias para armazenar as matrizes factores [L] (é inútil armazenar as matrizes [RJ], visto 
serem as transpostas de [L], é em geral inferior ao número de posições necessárias para armazenar a 
matriz [A] completa, podendo no entanto ser ligeiramente superior. Esta propriedade pressupõe a utiliza- 
ção de uma ordenação eficiente das colunas da matriz [AJ]. 

O algoritmo de cálculo digital decorre com a lógica seguinte. A partir do vector ILIM determina-se 


qual o primeiro nó a ser eliminado, isto é, quais as primeiras matrizes factores [1/2] e [rr] a Ser 
calculadas. Os elementos respectivos são obtidos pela equação (20) e armazenados nas posições previamente 
ocupadas pelos elementos da coluna correspondente âquele nó. 

Seguidamente são efectuadas todas as alterações nos elementos [B'] e [B'] devidos à eliminação do 
primeiro nó, de acordo com a equação (23). Os eventuais elementos novos são armazenados nas posições 
livres criadas durante a ordenação. O vector ILIM indica agora qual o segundo nó a ser eliminado, cal- 
culando-se [3] e [ pt] e repetindo-se a sequência até ao processamento do último nó. 

Esta lógica sofre apenas alteração devido à presença de nós tipo PU. Assim durante a bifactorização 
de [B"], quer os elementos das colunas correspondentes a nós tipo PU, quer os elementos das colunas tipo PQ 
com elementos numa linha tipo PU não são processados, 


A matriz [A], usada exemplo na secção 3.2, transforma-se, após a bifactorização em 


0,5 05 1 0 

0,5 04 0,4 1,8 
1 04 0384 | 2,54 
a e TS 
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São apenas armazenados os elementos não nulos (excluindo os elementos unitários do diagonal) das 
matrizes factores [L], situadas abaixo da linha a traço interrompido. As tabelas 1 e 2 da secção 3.2 trans- 
formam-se em; 


TABELA 1 

I 1 2 3 4 
DIAG(I) 05 04 0,384 — 0,058 
ILOC(I) 2 9 1 0 
NANU(I) 1 3 6 0 


ILIM (T) 


TABELA 2 
I 1 2 3 d o 6 f 8 9 10 
VALOR (K) 0,5 q 0,4 1,6 O 2,538 O ) 0 ) 
ILINHA (K) 2 3 3 4 2 2 3 ) O) O 
Ud) 4 O b o) 5] t 10 Ú 


IPROX (K) 2 


3.4 — Ordenação 


No programa elaborado utilizou-se um método da ordenação dinâmica da sequência de eliminação, com 
vista a minimizar a criação de novos elementos não nulos, 

O método consiste basicamente em escolher para pivot, em cada passo da eliminação de linhas e colunas 
de [B'] e [B'], usando o método de Gauss-Jordan, o elemento diagonal da coluna que tenha o menor 
número de elementos não nulos. Se mais de uma coluna estiver nestas condições, escolhe-se a coluna com 
número mais baixo. A eliminação de uma linha e uma coluna de idêntica ordem numa matriz (simétrica) 
é equivalente à eliminação do nó correspondente no grafo associado, pelo que se usarão indistintivamente 
as duas designações, Em relação à matriz [A] usada como exemplo, verifica-se por inspecção das tabelas 
da secção 3.2 que o primeiro elemento pivot é o elemento diagonal da coluna 1, visto que o vector NANU 
tem para todos os elementos o valor 2. 

A escolha do próximo elemento pivot depende do número de elementos não nulos em cada coluna 
da matriz a factorizar, após a eliminação da primeira linha e coluna, Devido à criação de novos elementos 
não nulos que esta operação em geral acarreta, há que simular a eliminação, sem contudo se calcular os valo- 
res dos novos elementos. Nesta fase, o que importa saber é se estes elementos existem ou não. 

Para determinar quais as colunas que vão ser reduzidas com a coluna a que pertence o elemento 
pivot examinam-se na tabela 2 os elementos não diagonais desta coluna. O índice da linha do primeiro 
elemento não nulo (dado pelo vector LINHA) indica qual a primeira coluna a ser reduzida (*). Conside- 
rando a matriz [A] como exemplo, esta seria a coluna 2, 


Para simular o processo de redução, comparam-se os índices das linhas dos elementos não diagonais 
da coluna a reduzir com a coluna redutora (pivot). Se numa dada linha o elemento da coluna redutora 
for não nulo e o elemento da coluna reduzida for nulo, há criação de um elemento novo. 

Após a simulação da redução da primeira coluna, procura-se o próximo elemento não diagonal da coluna 
redutora. O índice da linha deste elemento indica qual a próxima coluna a reduzir, O processo repete-se até 
que se esgotarem os elementos não nulos da coluna redutora. Em seguida determina-se por inspecção do 
vector NANU qual o próximo elemento diagonal que vai servir de pivot. Note-se que NANU é alterado 
durante o processo de simulação da ordenação devido à anulação dos elementos da linha a que pertence o 
pivot e à criação de novos elementos, 

No processo de eliminação de Gauss-Jordan, a linha e coluna a que pertence o elemento pivot são 
eliminadas e o elemento correspondente do vector independente modificado de acordo com as operações 


(*) Está-se apenas a considerar o caso de matrizes simétricas, 
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efectuadas. Na bifactorização o vector independente não é alterado, sendo as operações envolvidas na elimi- 
nação memorizados nas matrizes factores [R] e [L]. Os elementos destas matrizes vão ser armazenados 
nas posições previamente ocupadas pelas linhas e colunas eliminadas. Note-se que, se a matriz a factorizar 
for simétrica, apenas são armazenados os elementos das matrizes [LJ]. Durante a simulação da ordenação 
são anulados os elementos da linha a que corresponde o elemento pivot, de forma a criar posições livres na 
tabela 2 para armazenar os novos elementos não nulos. 


A presença dos nós tipo PU obriga a alterações na lógica que acaba de ser exposta, Como já foi 
referido na secção 3.2, os elementos das linhas e colunas correspondentes a estes nós não são alterados na 
simulação da matriz [B'"], sendo a sua posição nas tabelas dada pelos vectores LOC e NPROX, 

Para o cálculo dos elementos das matrizes factores [L] durante a factorização, os nós devem ser 
processados pela ordem determinada no processo de ordenação. É portanto necessário criar um novo vector, 
ILIM, na tabela 1, onde fica armazenada a ordem de eliminação determinada durante a ordenação. O índice 
de ILIM indica a ordem de eliminação e o seu valor indica o número do nó a ser eliminado: assim, se 
ILIM(n) = k, o nó k será o n-ésimo nó a ser eliminado. 


3.5 — Cálculo dos acréscimos [45] e [AU] 


As soluções [45] e [AU] das equações (8) e (9) obtém-se, de acordo com a equação (24), pela multi- 
plicação sequencial das matrizes factores de [B'] e [B"] pelos vectores [AP/U] e [4Q/U], respectivamente 


La [Re] Le] SI 1 Te E] [9] Lagço] (25) 
[4u| = [8º] [Rr] as [RD] [nt] EA Erad [1"2] [aQ/U | (26) 


Usando as técnicas de programação com esparsidade mencionados na secção 3.2, esta multiplicação é 


levada a cabo com extrema rapidez. Considere-se a multiplicação de um vector [8] pela matriz [299]: 
Tem-se 


(27) 


onde 


vo (28) 
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Do ponto de vista de cálculo a operação representada pela equação (27) decorre do seguinte modo. 

De acordo com ILIM (ver secção 3.4) determina-se qual o primeiro nó para o qual se vão calcular [45] e 
1 

[Au]. A partir do conhecimento deste nó procuram-se os elementos das primeiras matrizes factores ES 1 

e 48] alteram-se os valores de [AP/U] e [AQ/U] de acordo com a eg. (28). Em seguida determina-se 


através de ILIM qual o próximo nó a ser tratado e procuram-se os elementos de E raul e [1,027 alte- 
rando-se novamente os valores de [AP/U] e [A4Q/U]. O processo repete-se até ao último nó da rede. Salien- 
ta-se que as colunas das matrizes [L] são esparsas pelo que apenas alguns elementos dos vectores [4P/U] 
são alterados em cada multiplicação pelas matrizes [L'] e [L'] respectivamente, 


O produto de um vector [8º] por uma matriz factor ER] resulta num vector [BR] dado por 


t + | + mm = 


0 1 
po ) po j 
(K =1) (K=41) (K=1) 
j á mn 2 in : 
(K—1) (K=4) , (K—1) b (0) b (1) 
1 E 1 ê à (29) 
1 
O (1) 
po ) b, 
(0) (1) 
1 by, b, 
onde 
(1) (0) ' 
D, = b, La] 
n ” (K=1) 
! ( 
a (K-1) 
= a 


Já 


Dado que [R$] li [1 IE àparte o elemento diagonal, a multiplicação expressa pela equação (29) 
processa-se de um modo análogo ao descrito acima, com a diferença que apenas um elemento dos vectores 
[AP/U] e [AQ/U] é alterado em cada multiplicação, de acordo com a equação (30). Os elementos das matri- 
zes [R] são procuradas nas tabelas de armazenamento em sentido inverso ao da procura das matrizes [Lj], 
isto é, a partir do último nó a ser eliminado no processo de factorização. 


3.6 — Limites de potência reactiva nos nós tipo PU 


Depois de completada a segunda iteração, é feito um teste à potência reactiva calculada nos nós tipo 
PU, comparando-a com os limites máximo e mínimo especificados, Se um destes limites for violado a potên- 
cia reactiva no nó passa a ser constante e igual a este limite, isto é, o nó passa a ser do tipo PQ. Visto a 
convergência ser oscilatória, torna-se necessário prever a hipótese de a potência reactiva calculada voltar a 
entrar dentro dos limites numa iteração subsquente, passando neste caso o nó a ser novamente do tipo PU. 


O tratamento dos nós tipo PU em que os limites de potência reactiva são infrigidos pode ser feito por 
dois métodos: 


— O primeiro método consiste em acrescentar a [B"] a coluna e linha correspondentes ao nó PU em 
que a potência reactiva excedeu os limites e proceder a uma nova bifactorização. 
-— O segundo método consiste em corrigir a tensão por meio de um acréscimo AU, em cada uma das 


iterações subsquentes, de forma a anular o erro 4Q, = Q. — Qu O acréscimo d4U, relaciona-se 
com o erro AQ, por meio de um factor de sensibilidade S tal que 
AQ, 
AU =5 - 
= Mg U, (31) 
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S, é o elemento diagonal correspondente ao nó K da matriz das impedâncias nodais, podendo admi- 
tir-se aproximadamente que esta é a inversa da matriz [B"] aumentada com a linha e a coluna referentes 
ao nó K. S, é calculado apenas uma vez para os nós em que os limites foram violados, sem refactorizar a 
matriz [B"]. 

No programa desenvolvido, foi utilizado o segundo método, embora necessite de um número mais elevado 
de iterações, é, do ponto de vista de computação mais eficiente, já que não exige nova factorização de [B"]. 

Se num nó tipo PU os limites de potência reactiva são violados, esse nó é numerado em último lugar 
e a matriz [B”] aumentada será 


O a uia am 


Os elementos da coluna [C] e o elemento diagonal c | são detectados nas tabelas 2 e 1, respectiva- 
mente, por meio dos vectores LOC e NPROX, Estes elementos não foram alterados durante a ordenação 
e bifactorização. 

A equação (9) escreve-se nestas condições 


[4Q/U] | = | [BJ | IC] [4U] 


| | (33) 
AQ,/U, [Ji AU 


Se este sistema de equações for resolvido com recurso a uma eliminação gaussiana, o elemento dia- 


gonal c | que resulta de c | é o inverso do elemento correspondente da matriz inversa de [B"] ui . O ele- 


til 
+ 

mento c' | pode ser calculado facilmente sem inverter [B"] se se atender a que o produto das matrizes 

factores [L] de [B"] pela coluna [C] resulta numa coluna [C'] dada por: 


[e]= [1] [18-07 [109] [o] a 


que é idêntica à que resultaria de [C] por meio da eliminação gaussiana da matriz [B", Cl. Após esta 


eliminação a matriz [B"]” transforma-se em 


[B']" é [0] 
| e  spncerasr (35) 
[C 1 Can 
onde 
” 7 t 
[e] e [Sin Con  moE Cn: n ] (36) 


RE 
e [B"] é a matriz triangular superior que resulta de [B"]. Eliminando a última linha da matriz (35), o 


elemento RR transforma-se em 


n-1 2 
| E 
Can E “an = 1 (37) 
j=1 di 


= (38) 
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O erro àQ, é calculado no fim de cada iteração usando os valores das tensões em todos os nós cal- 
culados nessa iteração. A correcção AU, é calculada pela equação (31) e somada à tensão corrente nesse nó. 
Na primeira vez que é aplicada esta tensão é a especificada; nas iterações subsequentes será o valor corri- 
gido na iteração anterior. Se, numa dada iteração, a tensão calculada toma um valor superior à tensão 
especificada (se o limite de potência reactiva violado for o superior) ou inferior a essa tensão (se o limite 
violado for o inferior), a correcção deixa de ser implementada, passando o nó a ser tratado como um nó 
tipo PU normal. A potência reactiva calculada (usando a tensão especificada no nó K e nos restantes nós 
tipo PU em que os limites não foram violados) deve no entanto continuar a ser comparada com os limites, 
retomando-se o processo correctivo, caso estes sejam ultrapassados. 

Ainda que em dois ou mais nós sejam infrigidos os limites de potência reactiva, o factor de sensibili- 
dade é calculado para cada um deles sem atender a se existem ou não outros nós nessas condições, As linhas 
e coluna correspondentes a esse nó são por conseguinte acrescentados à matriz [B”"] base, 


3.7 — Transformadores com regulação automática em carga 


A partir da segunda iteração é também efectuado um teste à tensão nos barramentos com tensão 
controlada através de transformadores com mudança de tomadas automática em carga. 
O algoritmo de ajuste da relação de transformação consiste basicamente no seguinte. Se numa dada 


k 
iteração a relação de transformação do transformador que controla a tensão no barramento K for m e 


se a tensão calculada U, violar os limites máximo ou mínimo especificados, a relação de transformação é 


k 
corrigida do acréscimo Am, o dado pela expressão 


k 
am, O =t a(U-U,") (39) 


onde a é um factor empírico com as características idênticas às do ganho de uma cadeia de controlo com 
retroacção. No programa elaborado, um valor «=0,7 conduziu à mais rápida convergência, O valor corrigido 
da relação de transformação será 


K+1 K 


O sinal que precede w na equação (39) depende da localização do nó em que a tensão é controlada. 
Se este se situa do lado do enrolamento que possui as tomadas o sinal a utilizar deverá ser negativo, no 
caso contrário será positivo. 

Após a correcção de m,, é necessário alterar os elementos da matriz de admitâncias nodais, em virtude 
da modificação resultante nos esquemas equivalentes em « dos transformadores afectados. Em rigor também 
os elementos da matriz [B”] deveriam ser alterados, o que implicaria uma nova bifactorização. Verifica-se 
porém que, se a matriz [B”"] permanecer constante a velocidade de convergência do método é pouco afec- 
tada [10]. Dado que neste caso não é necessário refactorizar [B'"], a economia em tempo de cálculo é signi- 
ficativa. Os resultados obtidos pelos autores corroboram estas considerações. 

A variação da relação de transformação supõe-se contínua até à convergência total, após o que se fixa 
no valor real mais próximo do calculado e se procede à continuação do processo iterativo até nova conver- 
gência sem alteração da relação de transformação. Note-se que, com este procedimento, é possível que em 
alguns barramentos as tensões estejam fora dos limites especificados; em geral estes erros são desprezáveis. 


mun max 
Se m, em, forem respectivamente as relações de transformação mínima e máxima e Am, 
for o incremento, o valor real da relação de transformação mais próxima da calculada é 


min 
m =m, +K am, (41) 
onde K e o inteiro mais próximo de 
cal min 
m, = m, 
FK=——— (42) 
Am 


por defeito ou por excesso consoante | FK-K | <0,5 ou | FK-K | > 0,5 respectivamente. 
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4 — RESULTADOS NUMÉRICOS 


O programa de cálculo descrito na secção anterior foi utilizado na solução do fluxo de potência em 
diversas redes, nomeadamente as redes teste de 14 e 30 barramentos do IEEE, a rede teste de 56 barra- 
mentos da EDP (18.12.74 às 11.30 horas) e a rede de 63 kV da CRGE de 70 barramentos (situação de ponta 
em Janeiro de 1976). Na rede de 30 barramentos do IEEE consideraram-se três casos: sem violação dos 
limites de potência reactiva, com uma, e com duas violações. Na rede de 56 barramentos da EDP consi- 
deraram-se dois casos; Com e sem regulação automática em carga dos transformadores, Na rede de 70 
barramentos da CRGE, que inclui extensos percursos em cabo subterrâneo, estudou-se apenas o caso de trans- 
formadores com regulação automática em carga. 


Os cálculos foram executados no computador IBM 360/44 do Centro de Cálculo das Universidades 
de Lisboa, usando precisão simples. A tolerância admitida para os erros de fecho foi de 0,05 MWMVAR,; a 
redução desta tolerância para 0,01 MW/MVAR não revelou qualquer diferença nos valores das tensões (ex- 
pressas em p.u. com 4 casas decimais), 

Na tabela 1 apresentam-se alguns dados com interesse referentes às redes estudadas: valores máxi- 
mos de | AP | sede 8 | àQ | max após a primeira iteração e após a convergência, e o número de iterações 
necessário com o método do desacoplamento. Refere-se ainda o número de iterações usando o método de 
Newton-Raphson clássico, para efeitos de comparação. Dado que o processo iterativo com o método do desa- 
coplamento pode terminar após uma solução da equação (8) ou da equação (9), o número total de iterações 
pode ser inteiro ou fraccionário. 


TABELA 1 


Após iter. 1/1,5 


E die io Begui: Convergência vi de iterações 

Rede com Q lim autom. | | | Tai tói 

violado trans. [E | usas | JAO | mes | SP | mes | | AM | qraz Desacop. Faphada 
| TEEE [Ê | NÃO 47,24 T,713 | 0,0063 | 0,0018 
IEEE | | ) | NÃO 2,2 —asmo 0,0040 “0,061 
IEEE | 1 NÃO 4 oT7l “00014 O,QOTI 
| IEEE 2 | não | 2382 4371 | 00122 | 0,0436 
E DP Ê NÃO 513,47 5830 | 0,0454 0572 
E DP | | UH) | SIM o13,47 5830 0,0285 | 0,0229 
CRGE Ê Í do2,08 11,57 | O0178 0,0275 


NOTA: Valores de |A4P| e |4Q] em MW e MVAR, 


Do exame dos valores da tabela 1 pode concluir-se que, do ponto de vista do número de iterações neces- 
sário para atingir a convergência, o método do desacoplamento não se afasta sensivelmente do método de 
Newton-Raphson. Para a rede do IEEE de 30 barramentos com duas infracções de limites de potência reac- 
tiva, o número de iterações pelo método do desacoplamento é mesmo inferior ao número de iterações pelo 
método de Newton-Raphson, Dado que cada iteração pelo primeiro método é muito mais rápida que pelo 
segundo, o tempo total de cálculo é substancialmente menor pelo método do desacoplamento. 

Na tabela 2 apresenta-se uma decomposição do tempo total de cálculo, pelo método do desacoplamen 
to, para a rede do IEEE de 30 nós, sem violação dos limites de potência reactiva. 

De salientar os tempos extremamente curtos para a ordenação e bifactorização. Os tempos de leitura 
de dados e de formação das matrizes [Y], [B] e [B"] é comparativamente maior. Excluindo o tempo tomado 
pela saída de resultados que constitui 89 % do tempo total de duração do programa, pode verificar-se que 
o cálculo iterativo das tensões (3.5 iterações) demora apenas (0,6s, isto é, 15% do tempo total de cálculo. 
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TABELA 2 


Módulo 


— Leitura de dados e formação de [Y], [B] e [B'] 
— Ordenação 

—. Bifactorização 

— 3.5 iterações 


— Saída de resultados 


— Total 


5 — CONCLUSÕES 


O método do desacoplamento apresenta-se presentemente como uma alternativa superior ao método 
de Newton-Raphson, para a solução do fluxo de potência numa rede de energia, quer em exigências de memó- 
ria, quer em velocidade de cálculo, | 

Neste artigo dá-se conta do programa desenvolvido no Centro de Electrotecnia das Universidades de 
Lisboa utilizando este método, designadamente das técnicas de programação apropriadas às matrizes espar- 
sas, cuja aplicação ao cálculo de redes eléctricas permite reduções espectaculares da capacidade de memó- 
ria do computador e do tempo de cálculo. O uso da ordenação dinâmica das linhas e colunas das matrizes 
[B] e [B"] combinada com a bifactorização revelou-se extremamente eficiente do ponto de vista de veloci- 
dade de cálculo e simplicidade de programação quando comparado com outros métodos, designadamente a 
transformação [L] [D] [HJ], usada por Stott. 

Os resultados obtidos para as diversas redes testadas, embora ainda limitadas, revelaram-se satisfa- 
tórios, nomeadamente quanto à verificação da superioridade do método do desacoplamento em relação ao 
método de Newton-Raphson clássico. O número de iterações, para a mesma tolerância, é apenas ligeiramente 
superior na maioria dos casos, sendo igual num caso e mesmo inferior num outro, 

Presentemente os autores investigam a aplicação do método à análise de contingências, 
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7 — ANEXOS 


7.1 — Simbologia 


By, — elemento da matriz [AJ]. 

b, = elemento da matriz [B']. 

bi; — elemento da matriz [B']. 

B,, — parte imaginária de [X,;]. 

G,; — parte real de [Y,,). 

H, = elemento da submatriz [HJ]. 

J, j = elemento da submatriz [JJ]. 

L, E elemento da submatriz [L/|]. 

N e = elemento da submatriz |[NJ]. 

m, — relação de transformação do transformador i. 
n — número de barramentos da rede. 

P, — potência activa injectada no barramento i. 

Q, — potência reactiva injectada no barramento i. 

KR, Eme resistncia do ramo i—j. 

S. — factor de sensibilidade do barramento K. 

U, — módulo da tensão V.. 

V, — tensão no barramento i. 

a reactância do ramo i—j. 

Y,, — elemento da matriz de admitâncias nodais [TJ]. 
ô, — argumento da tensão V.. 

AS, — correcção do argumento da tensão 5,. 

AP, — erro de fecho da potência activa no barramento i. 
AQ, — erro de fecho da potência reactiva no barramento i. 
AU, — correcção do módulo da tensão U.. 


7.2 — Formulação do método do desacoplamento 
Retomem-se as equações (6) e (7) 
[AP] = [H] [ 43 1 (6) 


[49] = [L] [4U/U] (7) 
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onde os elementos de [H] e [L] são dados pelas equações (4a), (5a) e (9b). 
Atendendo a que, na generalidade das redes de energia, são válidas as seguintes aproximações: 


as equações (4a), (5a) e (Db) simplificam-se para 


H=L,=-U, B,, U, 1) (43) 


HEme=-U Bl, (44) 


li 


A linha i das equações (6) e (7) será 


n 
AP =- 5 U BjU, ds, (45) 
j=2 
n 
Q=- 5 UB,U,A4U, (46) 
i=2 = A ui td 
ou ainda 
aP, Ee 
——=- à BU, 48, (47) 
U. LE d d 
| = 
n 
AQ | 
—=- 3 B,4U, (48) 
U j Í 
! = 


NOTA: Admite-se que existe apenas um nó de balanço, que recebe o n.º 1. 
Sob forma matricial compacta as equações (45) e (46) escrevem-se 


[SP/U] = [B'] [UI] [48] (49) 
[50/U] = [B"] [4U] (50) 


Eliminando a influência do fluxo de potência reactiva no cálculo de AS (equação (49) ), para a que 
se igualam a 1.0 pu todos os elementos da matriz diagonal [U], as equaçes transformam-se finalmente em 


L4P/U] = [B'] [48] (51) 
[áQ/U] = [B"] [40] (52) 


Em princípio as matrizes [B'] e [B”] coincidem com a matriz [-B] (simétrica da matriz das sus- 
ceptâncias nodais [B]). Levando o processo do desacoplamento até as suas últimas consequências, não se 
consideram porém na formação de [B'] os elementos que afectam especialmente o trânsito de potência reac- 
tiva: capacidades ou induções em derivação (incluindo capacidades das linhas) e transformadores com relação 
de transformação variável, Também no cálculo de [B'] se despreza a resistência dos ramos da rede. Este 
ponto é de maior importância do que Stott sugere, pois a inclusão das resistências nos elementos de [B'] 
resulta num maior número de iterações até à convergência e pode mesmo conduzir à divergência do processo 
iterativo. 

Na formação da matriz [B”] ignora-se a eventual presença de desfasadores. Na matriz [B'] podem 
ser consideradas, o que, porém, resulta na destruição da simetria numérica desta matriz. É portanto pre- 
ferivel representá-los por meio de métodos alternativos, 


382 TÉCNICA 445/446 


NÚMERO 445/446 DEZEMBRO 1977/IANEIRO 1978 


ANO LII/LIHI 


VOLUME XXXIX 
C. D. U. 622.2:519.3 


Optimização do beneficiamento de sistemas 


de jazidas minerais 


[Apropriação racional dos recursos minerais ) 


RESUMO 


No planeamento a níveis nacional ou regional e 
também no âmbito de uma grande empresa responsá- 
vel por vários depósitos minerais, é-se conduzido ao 
problema da optimização global de sistemas de jaxzi- 
das, os quais podem conter minérios diferentes, 

Tomando como constrangimentos a taxa anual 
de produção de minério tal-qual e o teor minimo dos 
concentrados aceite pelos utilizadores, o objectivo con- 
siste em optimizar duas funcionais simultaneamente: 
a perda total em metal e a margem monetária global 
não actualizada. 

Para resolver este problema, usa-se a «função 
característica», já apresentada noutra oportunidade, 
mediante um processo de «homogenização» que se 
traduz na utilização de «teores reduzidos», 

Deste modo, constroi-se a função caracteristica 
global do sistema e pode obter-se facilmente a mar- 
gem monetária optimizante que conduz à perda mi- 
nima em metal compativel com os constrangimentos. 

O programa CARACT que optimiza as duas men- 
cionadas funcionais é sumarizado e algumas das suas 
virtualidades analisadas, 


1— A FUNÇÃO CARACTERÍSTICA OPTIMAL 


É sabido [5] que a escolha de critérios de rejeição 
e correlativa definição da massa minerável de qual- 
quer jazida só têm significado quando baseados em 
blocos com alguma representatividade tecnológica no 
que concerne à futura exploração, 

Por outro lado, critérios de rejeição baseados no 
teor tal-qual enfermam de aproximação grosseira [2], 
que consiste em admitir que a recuperação em me- 
tal daquele não depende. 

Assim, a construção do critério econômico a opti- 
mizar deve considerar (excepto em casos triviais) os 
resultados previsíveis da concentração e, se possível, 


JOSÉ QUINTINO ROGADO (1) 


ABSTRACT 


When planning for large areas, say a country, 
o mining district, or a great concern responsible for 
several mines, the exploitation and beneficiation poli- 
cies deal with complex systems, 

In a global view point, the best objective is to 
minimize the overall loss of metais and to maximize 
the total non-discounted monetary margin. 

The most important constraints framing the re- 
solution of such a problem are the minimum (or the 
maximum) grade levels of the concentrates to be 
produced and the annual rates of production admissi- 
ble by the market. 

The method developed in this paper is a quite 
general one and uses the concept of characteristic 
function elsewhere developed by means of an homo- 
genization process which leads to consider the de- 
posits set as q sole ore body, even if far apart ore 
bodies there exist. 

The computer program CARACT (built up for 
the a.m. method) is outlined and results already 
obtained with il are summarized, 

Finally, among other suggestions, the possibility 
of using this method for the evaluation of mineral 
resources in world scale, is advanced. 


da própria metalurgia ou processamento químico sub- 
sequente. Reconhece-se assim a necessidade de expli- 
citar sucessiva integração vertical na elaboração das 
políticas de exploração da maioria das jazidas, quan- 
to mais não seja para racionalizar as tradicionais 
«fórmulas de venda». 

A importância do bloco de exploração, acima in- 
dicada, conduz, já por si, à partição da jazida que, 
deste modo, se torna sistema discreto. E, quando 
os resultados da mineralurgia forem tomados em con- 
sideração, serão associações de blocos de exploração 
(constituindo blocos mineralúrgicos) que há que con- 


(1) Professor Catedrático do Instituto Superior Técnico. Departamento de Mineralurgia e Planeamento Mineiro. 
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siderar, já que a covariância dos parâmetros de con- 
trolo mineralúrgico é mais forte do que a existente 
entre os teores tal-qual, 

O problema que então se põe é encontrar para 
cada um daqueles blocos a forma mais adequada de 
beneficiamento. 

É cada vez melhor aceite que a política desejável 
consiste em maximizar a margem monetária global 
não actualizada, minimizando, simultaneamente, a 
perda em metal (ou em metais) de valor [5]. 

Se a reserva total Q «eventualmente» explorável 
for conhecida e se a taxa anual de extracção r for 
fixada, tal margem EB serã (v. lista de símbolos) 


B =Up—(fC + M(Q,r) + L(Q,r)) (1) 


Mostra-se [3], quando as alternativas de bene- 
ficiamento conduzem a custos quase constantes, que 
o teor T do concentrado global é atingido com perdas 
em metal U e monetária mínimas se a funcional 


*T 
I= [ | dT/dU | dT (2) 
A To 


for minimizada; o que corresponde a encontrar, con- 
siderando o conjunto 


A=[E:E= AR/A (R t) | (3) 


o sup E C A cuja configuração de tratamento dos 
painéis da jazida conduz ao valor de T fixado. 

Variando 'T, constroi-se no espaço UXC, a «curva 
característica (optimal)» da jazida na vizinhança de 
determinado processo de tratamento (ou de classe 
estreita (1) de processos de tratamento), 

E, se nos planeamentos sucessivos que se forem 
fazendo, os processos tecnológicos houverem sido bem 
escolhidos, em face do progresso técnico coevo e das 
características então conhecidas da jazida, as corres- 
pondentes curvas características podem gerar curvas 
de Lasky: Na época em que foram construídas tradu- 
ziam a solução optimal ideal para o aproveitamento 
da jazida, a qual não dependia dos valores de Q e r 
que se vieram a fixar. 

O interesse de qualquer projecto de beneficia- 
mento que se caracteriza por M(Q,r) + L(Q,r), que 
sº insere no mercado através de p, e toma em conta 
a localização da jazida, graças a f, pode avaliar-se 
à custa da curva caracteristica da jazida, 

Com efeito, sendo a margem monetária global 
não actualizada expressa em unidades do preço pro- 
-rata do metal de valor, dada por 


BSS . 


basta procurar o sup (U—G!) para determinar a 


MA A 
tríada (U,C, 7) optimal. A viabilidade económica do 


projecto resultará, portanto, de ser sup (U —G) > 0. 

Se tal não ocorrer e não se desejar explorar mina 
subvencionada, haverá que alterar o projecto ou re- 
jeitar blocos mineralúrgicos, 

A análise da política de rejeição conduz à conclu- 
são que a solução optimal (com minima perda em 
metal) consiste em abandonar, em prioridade, os blo- 
cos que conduzirão a menores teores de concentra- 
do [1]. 


2— O PROCESSO DE HOMOGENIZAÇÃO 


2.1 — Quando se ultrapassa a micro-análise, por 
exemplo, quando um sistema constituído por jazidas 
diferentes tem de ser considerado como um todo, a 
maximização da margem monetária global não actua- 
lizada e a simultânea minimização da perda em metal 
tornam-se problemas mais intrincados. Situações 
como esta ocorrem quando o projectista (ou o ges- 
tor) encara o planeamento a nível nacional ou regio- 
nal ou no âmbito de grande empresa concessionária 
de vários depósitos minerais, 

Hã, contudo, em geral, certa informação, na fase 
em que o projectista intervem, a qual permite, até 
certo ponto, enquadrar a situação. Assim, são conhe- 
cidas as reservas minerais e a taxa anual de produção 
futura de cada concentrado pode ser fixada com mar- 
gem de indeterminação não excessiva, depois de es- 
tudo individual das próprias jazidas, feita triagem 
de uns tantos processos de beneficiamento e após a 
análise dos mercados absorvedores dos futuros con- 
centrados, Deste modo, porque fica definido aproxi- 
madamente o rendimento em peso, pode fixar-se a 
taxa anual individual de mineração e, consequente- 
mente, a vida de cada depósito. 

Por outro lado, o projectista deve admitir que a 
política que vai adoptar será revista periodicamente 
(de 5 em 5 anos, de 10 em 10 anos, etc.) e ainda 
sempre que a conjuntura económica e/ou social hou- 
verem claramente mudado. 

Mesmo assim, uma das maiores dificuldades con- 
siste em comparar o valor de venda e, consequente- 
mente, o valor das perdas dos diversos metais a pro- 
duzir pelo sistema de jazidas. 

Com efeito, nem sempre meros critérios econó- 
micos podem ser adoptados, A exaustão relativa das 
reservas disponíveis no sistema, a degradação dos 
teores dos respectivos minérios «in situ» e as suas 
relações com o senário tecnológico, porventura a im- 
portância estratégica de alguns metais, etc. são fac- 
tores de influência «vaga» que o projectista deverá 
tomar em consideração, o que, aliás, pode ser feito 
atribuindo ponderadores diferentes aos valores dos 
diferentes metais, 

Por outro lado, em cada época em que se fizer 
o planeamento, o progresso técnico e os «hábitos» 
industriais fixam os teores mínimos das espécies de 
valor e os teores máximos das espécies prejudiciais 


(1) Estreita para que o pressuposto da constância de L seja válido. 
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dos concentrados a produzir. Trata-se de constrangi- 
mentos muito fortes que, mesmo ao nível de planea- 
mento nacional, não podem facilmente ser mitigados. 
Importa, por isso, que o método de planeamento os 
possa tomar explicitamente em consideração, 

2.2 — Suponha-se então que o sistema S, cujo be- 
neficiamento se pretende planear, é constituído por n 
depósitos D.: S = (D, o RR (4,0). A cada depósito 
estão alocadas a reserva em minério tal-qual Q, a 
taxa anual de mineração r,, o custo global de mine- 
ração e tratamento M (Q,r) + L (Q,r) e, também, 
o valor por unidade de metal p, (pro-rata), além dos 
custos de transporte f, por unidade de peso do con- 
centrado até FOB ou CIF. 

Considerando a união de todos os depósitos U D, 
obtém-se novo conjunto H = (Vi ie (ln), je 
(1,m,)) constituído pelos blocos mineralúrgicos, a 
cada um dos quais se alocaram alguns pontos (sejam 
Sa 6) da função de lavabilidade R(t), à custa de tes- 
tes normalizados representativos do(s) processo(s) 
escolhido(s). 

Adoptem-se os ponderadores a: ie (1,n) para 
os preços de venda dos metais a produzir e um dos 
metais e uma das minas como referencial, Defina-se: 


E 
U, = ta, .P;/P)-U, 

à = 1 

1 e (1,n) 
C, = (£/£).6, 


como o peso «reduzido» do metal (contido) no concen- 
trado a produzir e o peso «reduzido» deste, respecti- 
vamente. 

Consequentemente, o peso reduzido global de me- 
tal será 


U = D(v, ie (1,m) ) 


e o peso reduzido global do concentrado 


= “ 
C = Blc :i € (1,0) 
Obviamente, os custos globais de mineração e trata- 
mento serão, respectivamente 
“ 
M = E(m, ie (1,n)) 


| FR 2 (1, “le (1,n)) 


Assim fazendo, uma vez que se haja escolhido um 
par ( R,, t,,) para cada V, ; & margem monetária 
global da exploração de S até à exaustão (em face 
das reservas conhecidas) será 


MB =Up—(Ct +M+D) (4) 


relação que é formalmente análoga a (1) e que per- 
mite definir um teor generalizado global para o con- 
junto dos concentrados a produzir, a saber 
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T=U0=t%, (A Ta dl Gt) 


A relação (4) traduz, assim, o processo de homo- 
genização a que vai recorrer-se, para efeitos de cál- 
culo. 


3 — GENERALIZAÇÃO DA FUNÇÃO CARACTERÍSTICA 
OPTIMAL 


Sl — Como se referiu, a função caracteristica 
optimal de uma jazida traduz as configurações de 
tratamento que conduzem, para dado T, à menor per- 
da em metal. 


pe — 2 
Considere-se agora o espaço auxiliar U X "(CR 
onde é 


= “A “ s— 
T= (1/D,) pa U ja, p,/ dá &, = Et. 


Nesse espaço, a maximização da funcional 


garante que qualquer T é atingido com a perda mi- 
nima em metal, Portanto, considerando agora o con- 
junto 


Te | E:E- A R/ A ts) | 


encontrando-se o sup E CA cuja configuração de 
tratamento dos painéis do sistema conduz ao valor 
fixado de T, assegura-se tal teor com perda minima. 

E, como é óbvio, variando T, constrói-se a rela- 
ção U(T) optimal, isto é, correspondente às perdas 
mínimas em U. 

Esta relação induz no espaço UXC a função 
característica optimal U (C) da forma que se referiu 
em 2.1, emutatis mutandi>. Esta função U(C) diz 
respeito a conjunto de processos de tratamento (ou 
de classes estreitas de processos de tratamento) e 
permite também avaliar o interesse de determina- 


do conjunto de projectos aplicáveis às jazidas do 
sistema. Com efeito, basta, para tanto, procurar o 


32 —Se for sup (U-—G) >0, o sistema pode 
ser explorado com margem não negativa. 

Ficam assim, para cada bloco de tratamento Vas 
de cada jazida D,, definidos, em condições macro- 
-optimais, o rendimento em peso R,, o teor t,, e o 
rendimento em metal n, j do concentrado respectivo, 
o que permite reconstituir para cada D, a correspon- 


dente função característica (macro-optimal) U (C,) 
e a margem monetária global não actualizada U —G,. 

Pode então dar-se o caso de algumas jazidas 
apresentarem margens negativas. Mas, mesmo nesta 
situação, o processo de homogenização atrás desen- 
volvido define a melhor forma de explorar cada de- 
pósito mineral integrado no sistema, tomando em 
consideração os constrangimentos, 

As jazidas com margens positivas criaram um 
fundo de acumulação que permite explorar o sistema 
com margem não negativa e tal fundo foi constituído 
garantindo a perda minima global possivel em metal, 
em face do progresso técnico coevo e do conjunto de 
preços de venda e de transporte então praticados. 

Esta é, sem duvida, uma das grandes vantagens 
da análise macro-econômica apresentada. 

do — Por outro lado, é habitual existirem limites 
(inferiores) superiores para os teores dos diferentes 
concentrados em espécies de valor (prejudiciais), os 
quais são fixados rigidamente pelo mercado dos uti- 
lizadores, o que traduz, como se mencionou, determi- 
nados hábitos tecnológicos aparentemente sempre jus- 
tificáveis. 


Assim sendo, a construção de U(C) não é livre. 
Hã que limitar as configurações mineralúrgicas au- 
torizadas à restrição de conduzirem, para cada de- 
pósito D,, a teores acima (abaixo) dos fixados como 
soleira (teto) de rejeição comercial, 


4 — RESUMO DO PROGRAMA «CARACT» 


41 — O princípio estruturante expresso pelo ope- 


rador sup E (fixado T) torna possível e prático a 
construção da função característica optimal de um 
sistema de jazidas recorrendo a computador, por dis- 
pensar o cálculo combinatório. 

A fig. 1 resume o programa CARACT, que si- 


multaneamente calcula a característica U(C), a mar- 


gem monetária optimal sup(U-G) e U, CO, T, 
G, e os blocos rejeitados para cada uma das jazidas, 
satisfazendo a vários constrangimentos, a saber: 


(ij teores dos concentrados individuais acima 
(abaixo) das soleiras (tetos): 

(ii) margens monetárias individuais condiciona- 
das (opcional); 

(iii) margem moneteária global duplamente 
condicionada. 


E evidente que o programa é também utilizável 


para uma jazida única, porquanto a ROTINA DE 
HOMOGENIZAÇÃO, nesse caso não actuará. 


9— VIAS DE APROPRIAÇÃO RACIONAL DOS RECUR- 
SOS MINERAIS 


A escolha de Y muito grande e de X=0 permite 


obter U(C) descondicionada da economia de apro- 
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priação, isto é, a função característica optimal ape- 
nas condicionada pela conjuntura tecnológica de ju- 
zante, a qual é expressa pelas soleiras (tetos) de 
teores, 

Tornadas as soleiras nulas (os tetos muito gran- 
des) constrói-se a função característica descondicio- 
nada daquela conjutura, 

Realizando os dois descondicionamentos, obtém-se 
a função característica optimal específica da tecnolo- 
eia de beneficiamento coeva. 

Se, à escala mundial, se empreendesse a constru- 
ção de função caracteristica global, era por esta úl- 
tima versão que importava iniciar a construção de 
políticas de preços e de soleiras (e tetos) de teores 
que assegurassem a melhor apropriação possivel dos 
recursos minerais do Planeta em face do progresso 
técnico disponivel. 

Também, em escala mais reduzida, por exemplo, 
em países em via de lançar ou rever a política de 
aproveitamento dos seus recursos minerais, ou em 
associações de produtores de minérios, a utilização 
das diversas versões da função característica optimal 
é sem dúvida, de utilidade extrema. A dificuldade 
do seu emprego reside apenas na constituição de ban- 
co de dados que permita a criteriosa construção de 
tal função e a sua revisão periódica, 

Assim, a consideração das funções caracteristicas 
é básico para a estruturação do que se pode chamar 
os programadores nacionais de empreendimentos mi- 
neiros [4]. 

Em mais detalhe, a função característica permite 
orientar as negociações para a fixação dos preços dos 
metais e das tarifas de transporte [5] porque expli- 
cita as consequências globais referentes à totalidade 
das reservas minerais afectadas em termos de perda 
em metal e margem monetária não actualizada e 
induz a representação (classificação das jazidas) 
como afixos do espaço [valor do metal recuperável 
maximal, V x margem monetária maximal, B) CR. 

Esta classificação sugerirã a necessidade de in- 
crementar a investigação e o desenvolvimento e tam- 
bém de alterar os condicionantes a juzante, no sen- 
tido de aumentar as margens das jazidas com menor 
relação BM. Serão tais jazidas as que, em priori- 
dade, devem ser chamadas a contribuir, no futuro, 
no necessário processo de acumulação. O método des- 
crito tem, portanto, virtualidades diacrônicas que 
poderão vir a orientar opções racionais na apropria- 
ção dos recursos minerais, 

No pormenor, refira-se apenas que o programa 
CARACT permite tratar rapidamente sistemas para 
os quais o uso de combinatória tornaria a utilização 
do ordenador extremamente dispendiosa ou impossí- 
vel. Num exemplo apresentado [5] encontra-se em 
15 segundos a função característica que, de outro 
modo calculada, exigiria que se lidasse com combi- 
natória de 37 820 configurações, 
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DADOS 


ROTINA DE HOMOGENIZAÇÃO 


CONSTRUIR E | 


CONSTRUIR A PRÓXIMA CONFIGU- 


RAÇÃO COM SUP E E CALCULAR G 


Mae ULTIMA 
CALCULAR U; C4 Ts; Gs 


CONFIGU- 
RAÇÃO 


- NÃO - SIM 


DAS JAZIDAS INDIVIDUAIS 


ULTIMA 


CONFIGU- NÃO 


RAÇÃO 
SIM 


CALCULAR SUP [U — 6] 


| ROTINA DE REJEI- 
ÇÃO DE BLOCOS 


IMPRIMIR U(G),SUP (U — G) Us Cpo Tio 


U.(C.), U.— G. E BLOCOS REGEITADOS 
L 1 L 1 


Fig. 1 
Fluxograma do programa CARACT 
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6 — LISTA DE SÍMBOLOS NÃO DEFINIDOS NO TEXTO 


6.1 — Referentes a uma única jazida 


— despesa total de mineração 


Pemtmo 


— despesa total de tratamento 


ij 


— taxa anual de mineração do tal-qual 
— teor do concentrado global 


— metal contido no concentrado global 


— peso do concentrado giobal 


o e MG 


em unidades «pro-rata» de metal 


6.2 — Referentes a sistemas de jazidas 


[1] 


E 
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— valor da unidade de metal («pro-rata» ) 


— frete por unidade de peso de concentrado 


— reserva total «eventualmente» explorável 


— margem monetária total não actualizada 


— rendimento em peso genérico; 7 rendimento em metal genérico 


— teor genérico do concentrado; A operador transferência [2] 


= 1/p. (no + M(Q, r) + L(Q, r)) despesas totais de transporte, lavra e tratamento expressas 


P, — valor da unidade do metal i 

Ê, — frete por unidade de peso do concentrado i 

P, — «pro-rata» de referência 

E — frete de referência 

U, — metal contido no concentrado global i 

e, — peso do concentrado global i; T, = U,/C, teor idem 

KR, — rendimento em peso do concentrado produzido pelo bloco mineralúrgico j da jazida i; t, teor 


do concentrado idem; n,, 
| 


rendimento em metal idem 


Es — teor do concentrado genérico <> t, | à, P;/P, 5/15 
w 


— denominação e massa do bloco mineralúrgico j da jazida i 


B, — margem monetária total não actualizada da jazida i 
G =1/p. (£f C+M + L) despesas totais de transporte, lavra a tratamento do sistema, expressas 


em unidades «pro-rata» do metal de referência 


X e Y — constantes arbitrárias positivas 
G, = 1/p, (£€, +M(Q,.r,) +I(Q,r,)) 


m — número de pontos de funcionamento mineralúrgico do bloco i 
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Resumo dos Artigos publicados na «Técnica» n.º: 445/446 
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C. D. U. 378.962 


E. R. DE ARANTES E OLIVEIRA 


Técnica N.º 445/4465 — LII/LII] — 12/1-1977/78, p. 129 a 138 


Consideram-se vários pontos básicos dos quais depende 
o futuro do Instituto Superior Técnico: os problemas da orga- 
nização geral do ensino da engenharia no Pais, do tipo e 
quantidade dos discentes, do corpo docente, da formação 
de professores através dos cursos de pós-graduação, do 
serviço docente em tempo inteiro, da organização da inves- 
tigação dentro da Escola e das relações com a Indústria. 


C. D. U. 378.962:62.001 
L. VALADARES TAVARES 


ENSINO DE ENGENHARIA E TRANSFORMAÇÃO 
DA SOCIEDADE 


Técnica N.º 445/446 — LI[LIHO — 121-1977/78, p. 139 a 142 


Neste trabalho analisa-se a acção profissional dos enge- 
nheiros no contexto dos processos de transformação da so- 
ciedade mostrando-se não ser actualmente possivel nem dese- 

— jável deduzir as suas necessidades de formação a partir de 
| um quadro estático de previsões feitas a curto ou a médio 
| prazo. Em contrapartida, traçam-se algumas linhas de evolução 
| julgadas mais importantes e a partir desta análise de tendên- 


QUE FUTURO PARA O 1.5.7.7? 
| 
| 
| 


nharia em Portugal. 
| Em conformidade com estas orientações apresentam-se 
| ainda algumas sugestões relativas ao conteúdo e à estrutura 
do ensino de Engenharia. 


C. D. U. 378.962 


e 


JOSE MANUEL FONSECA DE MOURA 


ENSINO E ESCOLAS DE ENGENHARIA: TÓPICOS 
PARA DISCUSSÃO 


Técnica N.º 445/4460 — LIILIH — 121-4977/78, p. 143 a 152 


Analisa-se o ensino e os tipos de escolas de engenharia, 
abordando quatro tópicos, nomeadamente a estrutura curri- 
cular, os graus e titulos outorgados, algumas questões peda- 
gógicas, as funções e o enquadramento do pessoal docente. 
Como principios orientadores têm-se a profissionalização da 
| escola, o carácter formativo do seu ensino e sua integração 
| no contexto social. Como metodologia distinguem-se na dis- 
cussão os niveis global e local. O primeiro corresponde aos 
objectivos básicos. O segundo permite a flexivel adaptação 
destes ao contexto. 
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cias apontam-se diversos objectivos para o ensino de Enge- 


C. D. U. 378.962:371.214 
JOSÉ MANUEL FONSECA DE MOURA 
PLANIFICAÇÃO DO CALENDÁRIO ESCOLAR 
Técnica N.º 445/4486 — LILI — 1214977078, p. 153 a 156 


Nesta comunicação estuda-se a organização do calendário 
do ano escolar das escolas de engenharia, considerando dois 
niveis de planificação. No global estabelece-se a divisão do 
ano escolar em dois grandes periodos — o ano lectivo e o 
periodo complementar. No local defende-se a quantificação 
da actividade escolar e, a partir desta, uma nova distribuição 
dos tempos lectivos, 


C D. U. 371.3:378.562 


MARIA DE LOURDES TEDESCHI| BETTENCOURT 
ALFREDO MANUEL FLORENÇO TOVAR DE LEMOS 


EXPERIÊNCIA EM PEDAGOGIA DE GRUPO 
Técnica N.º 445/446 — LILI — 121497778. p. 157 a 164 


Dá-se conta dos principios pedagógicos que constituem 
o apoio científico da pedagogia de grupo e caracteriza-se a 
situação circunstancial em que decorreu a experiência nos 
anos de 1970/71 e 72, envolvendo um total de cerca de 100 
alunos. Transcreve-se, em resumo, e a partir dos relatórios 
feitos pelos alunos participantes na experiência, as conclu- 
sões que defendem a pedagogia de grupo. como um processo 
natural de aprendizagem humana. Considerou-se como funda- 
mental nesta experiência a formação que permite a inserção 
dos futuros engenheiros na actividade profissional, onde se 
afirmem como técnicos de elevado nivel e como pessoas em 
perfeito e harmonioso equilibrio. 


C. D. U. 62.001.5 


PEDRO LOURTIE 
VENTURA SAMPAIO 
NELSON ESTEVES 
MARIO CORDEIRO 


DA INVESTIGAÇÃO AO PROJECTO DE ENGENHARIA 
Técnica N.º 445/446 — LI[ÍLIH — 12/1-1977/78. p. 165 a 168 
Faz-se uma análise da actual estrutura da Engenharia em 


Portugal, de investigação, desenvolvimento e projecto, apon- 
tando-se vias de solução de alguns dos defeitos encontrados. 


Synopsis of articles published in «Técnica» n.º: 445/446 


U, D. C. 378.962:371.214 
JOSÉ MANUEL FONSECA DE MOURA 
PLANNING OF THE SCHOOL YEAR CALLENDAR 
Técnica Nº 445446 — LIIILHI — 121-1977/78, p. 159 a 156 


The paper studies the engineering schools year callendar. 
It considers two levels of plannification. The global level 
establishes its division in two main periods — the school! 
year itself and a complementary period. The local level quan- 
tifies the achool related activity and from this, it advocates 
a new different organization of the present students schedule. 


U. D. C. 371.3.378.062 


MARIA DE LOURDES TEDESCHI BETTENCOURT 
ALFREDO MANUEL FLORENÇO TOVAR DE LEMOS 


EXPERIENCE IN TEAM TEACHING 
Técnica N.' 445446 — LI[LH| — 121497778, p. 157 a 164 


The pedagogical principles constituting the scientific 
support for group pedagogy are presented, and the circumstan- 
ces under which the experience, involving approximately 100 
pupils, was held in 1970, 1971, and 1972, are characterized. 
The conclusions, taken from the reports of the participating 
pupils, in defense of team teaching as a natural human lear- 
ning process, are summarized. A training process enabling 
future engineers to embark upon their professional careers 
both as highly qualifisd technicians and as well-balanced hu- 
man beings has bsen considered a fundamental factor in this 
experience 


U. D. C. 62.001.5 
PEDRO LOURTIE 
VENTURA SAMPAIO 
NELSON ESTEVES 
MARIO CORDEIRO 
FROM RESEARCH TO INDUSTRIAL DEVELOPMENT 
Técnica N.º 445/4456 — LIVLIH — 12/1-4977/78. p. 165 a 168 
The activities of the institutions In Portugal dedicated to 


research and development in the engineering fields Is briefly 
reviewed and a new plan of organisation is suggested. 


LIL / LIlI— December 1977 | January 1978 


U. D. C. 378.962 
E. R. DE ARANTES E OLIVEIRA 
WHICH FUTURE FOR |.8S.T.? 
Técnica N.º 445/4460 — LILI — 12/1-1977/78, p. 129 a 138 


The paper considers several basic points on which the 
future of IST will depend: the structure of engineering educa- 
tion in Portugal, the kind and number of engineering students, 
the qualitative and quantitative composition of the teaching 
staff, the organization of post-graduation courses, the full-and 
part-time teaching service, the organization of research, the 
relations with Industry. 


U. D. C. 378.962:62.001 
L. VALADARES TAVARES 


EDUCATION IM ENGINEERING RELATED TO SOCIETAL 
DYNAMICS 


Técnica N.º 445446 — LI/LHI — 12/1-1977/78, p. 139 a 142 


The impact of the engineers professional activity on the 
developing processes of portuguese society is studied and 
several remarks on the recommended methodology to set 
up the objectives for Education in Engineering are deduced 
from this analysis, 

Somme targets and suggestions about the structure and 
contents of engineering education in Portugal are also pre- 
sented in this paper. 


U. D. C. 378.962 
JOSE MANUEL FONSECA DE MOURA 


THE ENGINEERING SCHOOLS AND ITS TEACHING/LEAR- 
NING PROCESS: TOPICS FOR DISCUSSION 


Técnica N.º 445/446 — LII/LHI — 121-1977/78. p. 143 a 152 


The paper considers the teaching/learning process and the 
types of engineering schools. It analyses four topics, namely: 
1) the curricula structure, ll) the degrees and academic titles, 
Hi) the pedagogy, Iv) the functions and structure of the faculty. 
The underlying principles of the analysis are the professiona- 
lization of the schools, the formative caracter of ita teaching! 
learning process, and its integration in the social structure. 
As methodology two levels are distinguished: the global level 
and the local level. The first corresponds to the basical 
objectives. The second allows for flexible concretization of 
these objectives. 
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Resumo dos Artigos publicados na «Técnica» n.º: 445/446 


Ano LII / LIH — Dezembro 1977 | Janeiro 1978 


Tn a mma Es a e is O ss — — — = — O = e 


J. S. BRAZÃO FARINHA 
O ENSINO DA ENGENHARIA CIVIL NO NOSSO PAIS 
Técnica N.º 445/446 — LII/LII — 12/1-1977/78, p. 169 a 171 


Refere-se a posição actual do ensino da engenharia e 
indica-se uma necessidade de evolução, no sentido de poder- 
mos enfrentar os problemas concretos da comunidade nacio- 
nal, 

A primeira etapa da evolução a encarar é a preparação 
para a cobertura técnica do Pais, com vista aos problemas 
do quotidiano, dando satisfação às necessidades de resolução 
das questões básicas do nosso Povo, neste sector. 

Não sendo uma faceta única, ela é sem dúvida a etapa 
mais urgente e prioritária, abrangendo um maior número de 
técnicos a formar anualmente. 


C. D. U. 371.214:624 
ANTONIO DE CARVALHO QUINTELA 


C. D, U. 31.214:624 
PRINCIPIOS PARA A REESTRUTURAÇÃO DO ENSINO 
DA ENGENHARIA CIVIL 


| 
Técnica N.º 445/446 — LII/LII — 121-1977/78, p. 173 a 179 
| 


Caracterizam-se os perfis dos engenheiros civis e defi- 
| nem-se os principios que se entende deverem orientar a 
formação a facultar nos bacharelatos e nas licenciaturas, em 
| Engenharia Civil, e as respectivas condições de acesso. Tra- 
ta-se também do ensino de pós-graduação e faz-se uma 
breve referência à instituição de cursos para formar «técnicos 

| de Engenharias». 
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C. D. U. 371.214:622 
JOSE QUINTINO ROGADO 


NOTA SOBRE OS CURSOS UNIVERSITARIOS PORTUGUESES 
DE ENGENHARIA DE MINAS 


Técnica N.º 445/446 — LII/LHI — 121-A977/78, p. 181 a 186 


Analisa-se a natureza dos cursos universitários de Enge- 
nharia e detalha-se a situação dos de Mines em Portugal, 
sugerindo-se certos ajustamentos. 


CONSIDERAÇÕES SOBRE O TEMA: A FORMAÇÃO 
DO ENGENHEIRO 


Técnica N.º 445/4465 — LIWLII — 121-1977/78. p. 187 a 190 


Partindo dos conceitos de Formação e de Engenharia 
procura-se explicitar o que o autor pensa estar neles impli- 
cito, organizando, assim, a defesa de uma formação técnico- 
-sconômica e social para o engenheiro. ao longo da sua 
graduação. 

Propõe-se a abolição do «divórcios escola engenheiro atra- 
vês de cursos de pós-graduação para actualização e aprofun- 
damento de todos aqueles aspectos propostos na graduação. 

Faz-se uma referência muito especial à Informática, ten- 
tando diversas perspectivas para os beneficios que resultam 
do seu conhecimento e utilização. 

Finalmente faz-se referência ao caso de uma conceituada 
escola inglesa cuja orientação e actuação estão na linha 


C. D. U. 37.015:378.962 
H. SIMÕES SANTOS 
das propostas do autor, 


fi 


Es 
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U. D. G. 371.214:622 
JOSE QUINTINO ROGADO 


THE MINING ENGINEERING COURSES OF PORTUGUESE 
UNIVERSTIES 


Técnica N.º 445/446 — LIVIO — 12/1-4977/78, p. 181 a 186 


The nature of University Courses for Engineers is com- 
mented with Increase detail for the Portuguese Mining Cour- 
308. 

Consequently, some eadjustments in these courses are 
suggested 


U. D, C. 37.015:378.962 
H. SIMÕES SANTOS 
SOME HINTS ON ENGINEERING GLOBAL EDUCATION 
Técnica Nº 445446 — LILI — 121197778, p. 187 a 190 


The author attempts to point out what he considers to 
be implicit in the concepts of global education and engine- 
ring 

The abolition of the engincer's «divorce= from his college 
through post-graduate courses to up-dating and deepening of 
all those matters included in his graduation, is proposed. 

A special reference is made to informatics, attempting 
various perspectives resulting from its Knowledge and use. 

Finally, a reference is made to a renowed English college 
whose orientation and functioning are similar to the author's 
proposals. 
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U. D. C. 371.214:624 
|. 8. BRAZÃO FARINHA 
TEACHING OF CiVIL ENGINEERING IN PORTUGAL 
Técnica N.º 445/4465 — LIVLHI — 12/1-1997/78, p. 160 a 171 


it is refered the actual state of the engineering teaching 
and of an evolution necesaity, in order to enable to face the 
real problems of the national community. 

The first stage of the evolution to face is the preparation 
for the National technical cover responding to the necessities 
of solving our People fundamental questions, in this field, 
taking in consideration the daily problems. 

Although It is not an unique facet, it is, undoubtedly, the 
most urgent and with more priority stage. dealing with a 
larger number of technicians to graduate each year. 


U. D.C. 371.214:624 
ANTÔNIO DE CARVALHO QUINTELA 


PRINCIPLES FOR THE RE-ORGANIZATION OF CIVIL 
ENGINEERING EDUCATION 


Técnica N.º 445446 — LIH/LHI — 12/1-1977/78, p. 173 a 179 


The civil engincers profiles are characterized and the 
principles that should guide the schoolar curricula and admis- 
sion conditions are dicussed. Postgraduate studies are also 
considered and a brief reference to courses to prepare «techni- 
cianss is made. 
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C. D. U. 72.011.1 
RUY JOSE GOMES 


NECESSIDADES HUMANAS E EXIGÊNCIAS FUNCIONAIS 
DA HABITAÇÃO 


Técnica N.º 445/4846 — LIÍLIHIE — 120497778, p. 191 a 219 


Apresentem-se comentários de generalidade relativamente 
à definição das necessidades humanas e consequentes exi- 
gências funcionais da habitação e à sum confrontação com os 
condicionalismos de conjuntura que, na solução das carén- 
cias habitacionais, possam constranger ou adiar a satisfação 
plena daquelas exigências. 

Agrupam-se, por identidade de objectivos, e hlerarqui- 
zem-se, por grau de essencialidade, as exigências a consi- 
derar. 

Definem-se os principios exigenciais e, com apoio na 
documentação referenciada, caracterizam-se as exigências se- 
guintes: 

— Exigências de SEGURANÇA 

— ESTRUTURAL E CONSTRUTIVA 
— peronte solicitações correntes, de serviço 
— perante solicitações de ocorrência excepcional 
— perante solicitações de ocorrência acidental. 
em uso 
— contra RISCOS DE INCÊNDIO 
— na OCUPAÇÃO E USO 
— da circulação na habitação e acessos 
— contra quedas de andares elevados 
— das instalações e equipamentos 
— Exigências de SAUDE 
— por EQUILIBRIO TÉRMICO 
— por CONDIÇÕES DE SALUBRIDADE 
— de abastecimento de água 
— de evacuação de águas e esgotor 
— de higiene pessoal 
— de higiene doméstica 
— de pureza do ar embiente 
— de limpeza e desinfecção da habitação 
— Exigências de CONFORTO 
— por CONDIÇÕES ACUSTICAS 
— por CONDIÇÕES VISUAIS 
— por DISPONIBILIDADE DE ESPAÇO 
— Exigências de SATISFAÇÃO 
— por SEGURANCA ADICIONAL 
— na OCUPAÇÃO DOS ESPAÇOS 
— por CONDIÇÕES DE PRIVACIDADE 
— por FACILIDADES DE USO 
— por SUPLEMENTO DE CONFORTO 
— Exigências de ECONOMIA 
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C. D. U. 624.057 
JOSE A. TEIXEIRA TRIGO 


INDUSTRIALIZAÇÃO E QUALIDADE DA CONSTRUÇÃO 
Técnica N.º 445/446 — LII/LIIH — 12/1-1977/78, p. 221 a 230 


Analisa-se a forma como se vem processando a industria- 
lização da construção de edificios e as suas consequências 
sobre a qualidade da construção, com especial ênfase sobre 
a actual situação da construção em Portugal. Esta análise 
baseia-se na experiência pessoal adquirida no LNEC em con- 
tacto com numerosas obras e sistemas de construção indus- 
trializada, e nos estudos que, sobre estes temas, vem reali- 
zando no âmbito do Plano de Edifícios desenvolvido no mesmo 
Laboratório. 

Depois de se sintetizar a metodologia tradicional utilizada 
no controle de qualidade da construção, analisam-se as per- 
turbações nela introduzidas pelos novos processos de cons- 
trução industrializada. Referem-se as principais deficiências 
que se vêm encontrando. Propõem-se acções de apolo a 
novas metodologias de controle de qualidade, adequadas a 
uma indústria de construção renovada. 
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C. D. U. 699.81 


ALFREDO MANUEL FLORENÇO TOVAR DE LEMOS 
ILDEFONSO CABRITA NEVES 
JOAQUIM DA CONCEIÇÃO VALENTE 


ACÇÃO DOS INCÊNDIOS SOBRE A CONSTRUÇÃO 
Técnica N.º 445/4465 — LIHLIHI — 12/-1977/78, p. 231 a 238 


Chama-se a atenção para a necessidade de serem actuali- 
zados entre nós os conhecimentos relativos à acção dos 
incêndios sobre a construção, cujos aspectos principais se 
procurou pôr em evidência, no sentido de se obter um nivel 
de protecção adequado. 

Trabalha-se com os conceitos de resistência ao fogo dos 
elementos de construção relativamente às funções de com- 
partimentação e de estabilidade. Referem-se por fim a alguns 
aspectos econômicos ligados à problemática dos incêndios. 
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U. D. C. 699.81 


ALFREDO MANUEL FLORENÇO TOVAR DE LEMOS 
ILDEFONSO CABRITA NEVES 
JOAQUIM DA CONCEIÇÃO VALENTE 


FIRE AND CONSTRUCTION 
Técnica N.º 445/446 — LILI O — 12/1-1977/78, p. 231 a 258 


Atention is called to the need of bringing up to date in 
this country the knowledge about the action of fire on cons- 
truction, whose essential aspects are emphasized in order 
that an appropriate level of protection may be achieved. 

The concepts of fire resistance in compartment and stabi- 
lity functions are used. Finally, are refered some economical 
aspects involved In fire. 
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U. BD, O. 720141 
AUY JOSE GOMES 


HUMAN HEEDS AND FUNCTIONAL REQUIREMENTS 
FOR HOUSING 


Técnica N.º 445/4846 — LILI — 12/1-1977/78, p. 191 a 219 


Some general comments are presented on the definition 
of human needs and consequent functional requirements for 
housing, and on the confrontation with conjunctural limitations, 
which in the solution of housing needs may hamper or postpo- 
ne full satisfaction of these needs. 

The requirements to be considered are grouped according 
to the identity of their aims and graded according to their 
degree of priority. 

The requirement principles are defined and with the help 
of the documentation referred to the following requirements 
are characterized; 

— SAFETY requirements 

— STRUCTURAL AND CONSTRUCTIVE 
— under current service loads 
— under exceptional loads 
— under accidontal lozds, in use 
— against FIRE RISKS 
— under OCCUPATION and USE 
— of circulation in the house end its approaches 
— against falls from high sicreys 
— of installations and equipments 
— HEALTH requirements 
— through THERMAL EQUILIBRIUM 
— through SALUBRITY CONDITIONS 
— of water supply 
— of water and sewage disposal 
— of personal hygiene 
— of domestic hygiene 
— pf ambient air freshness 
— of cleanliness and disinfection of the house 
— COMFORT requirements 
— through ACOUSTIC CONDITIONS 
— through VISUAL CONDITIONS 
— through SPACE AVAILABILITY 
— Requirements as to SATISFACTION 
— through INCREASED SAFETY 
— in SPACE OCCUPATION 
— through PRIVACY CONDITIONS 
— through FACILITY OF USE 
— through ADDITIONAL COMFORT 
— ECONOMY requirements 


U. D. C. 624.057 
JOSE A. TEIXEIRA TRIGO 


BUILDING INDUSTRIALIZATION AND QUALITY 
Técnica N.º 445/4846 — LIIÍLII — 12/1-4977/78, p. 221 a 230 


Theo way in which the industrialization of building cons- 
truction ts evolving and how this is affecting the quality of 
construction are analysed, with particular emphasis for the 
present situation of construction in Portugal. 

This analysis is based on the authors experience at the 
LNEC. in contact with numerous works and industrialized 
construction systems, ând on the studies he has been carrying 
out within the framework of the Plan for Building Construction 
developed et this Laboratory. 

After synthesizing the methods usually adopted in the 
quality control of construction, the author analyses the changes 
introduced by the new methods of industrialized construction. 
The most important deficiencies are refered to. Actions for 
the support of new quality control methods are suggested, 
in view of a new and better industry of construction. 
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C. D. U. 624.07/012.4]8 
o. POMPEU DOS SANTOS 


LIGAÇÕES EM ESTRUTURAS PREFABRICADAS DE BETÃO 
Técnica N.º 445/446 — LIJ/LIM — 121-4977/78, p. 239 a 254 


à prefabricação das estruturas dos edificios implica a 
necessidade de executar em obra as ligações entre os ele- 
mentos prefabricados. À concepção, projecto e execução das 
ligações envolve a resolução de certas questões especificas 
normalmente não contempladas na regulamentação aplicável 
às estruturas de betão armado e (ou) pré-esforçado ou à 
estrutura metálica. 

Neste documento abordam-se as ligações entre as peças 
prefabricadas das estruturas reticuladas de betão, sob o ponto 
de vista das suas exigências, nomeadamente de segurança 
estrutural, de durabilidade, de montagem, etc., procurando-se 
apresentar os modos de como essas exigências poderão ser 
satisfeitas. Apresenta-se aínda uma classificação de ligações 
de acordo com vários parâmetros, fazem-se algumas conside- 
rações sobre a concepção e execução de alguns tipos de 
ligações, apresentam-se alguns exemplos de soluções realiza- 
das por processos diversos e finalmente faz-se a apresentação 
dum programa de estudos a desenvolver no âmbito deste 
tema. 


C. D. U. 627.8.09(679.4) 
ALBERTO ABECASIS MANZANARES 
ANTONIO CARRISSO 
ANTONIO DE CARVALHO QUINTELA 
ANTÔNIO JOSE DA COSTA SANTOS 
JOÃO DE SALVADOR FERNANDES 


A ENGENHARIA NACIONAL E À INSTALAÇÃO DE UNIDADES 
DE GRANDE POTÊNCIA. O CASO DE CABORA-BASSA 


Técnica N.º 445/446 — LI/LIN — 12/1-4977/78, p. 255 a 266 


Apresenta-se o caso de Cabora-Bassa como testemunho 
da capacidade da engenharia portuguesa para & concepção e O 
projecto da instalação de unidades hidroeléctricas de grande 
potência. 

Depois de uma descrição geral do aproveitamento, refe- 
rem-se resumidamente os principais condicionamentos e pro- 
blemas enfrentados, evidenciando a necessidade da sua con- 
sideração conjunta pelos fornecedores dos equipamentos e 
pelos projectistas das instalações. 


C. D. U. 621.313.12:519.28 


M. FERREIRA DE OLIVEIRA 
1. REIS CARDOSO 


ESTADOS DE AVARIA PARCIAL DE GRUPOS EM ESTUDOS 
DE RESERVA GIRANTE 


Técnica N.º 445/446 — LI/LH| — 121497778, p. 267 a 274 


A fixação da reserva girante feita em termos probabilis- 
ticos tem por objectivo garantir que a probabilidade da carga 
exceder a potência disponivel não é superior a um valor 
fixado como admissivel. 

Estudos desta natureza fazem-se conhecendo as probabi- 
lidades, como função do tempo, dos grupos do sistema se 
encontrarem nes «estados» em que podem residir; a cada 
um dos possiveis estados em que um grupo pode residir está 
associada uma potência máxima que o grupo é capaz de 
produzir. 

O presente trabalho apresenta uma metodologia que per- 
mite calcular as referidas probabilidades visando a sua utili- 
zação em estudos de fixação da reserva girante de um sistema 
electroprodutor. 


TÉCNICA 445/446 


Ano LH /LIH — Dezembro 1977 | Janeiro 1978 


C. D. U. 621.39(469) 


MANUEL JOSE ABREU FARO 
PEDRO LUIS BORGES TEIXEIRA 


CONCEPÇÃO E PROJECTO DE SISTEMAS 
DE TELECOMUNICAÇÕES NO PAIS 


Técnica N.º 445/446 — LII/LIH — 12/1-1977/78, p. 275 a 295 


C. D. U. 665.51/58 
LUIS L. RAMALHETE 
PETROQUIMICA 


Técnica N.º 445/446 — LII/ÍLII — 121497778, p. 297 a 303 


C. D. U. 665.51/.58 
JOÃO PEDRO ANTUNES GOMES 


O COMPLEXO PETROQUIMICO DE SINES; CONCEPÇÃO 
E ESTADO ACTUAL DOS TRABALHOS 


Técnica N.º 445/446 — LII/LIIH — 12/11-1977/78, p. 305 a 312 


Após se dar conta dos diplomas legais que lançaram 
a indústria petroquimica de olefinas, apresenta-se a concepçã! 
do Complexo Petroquímico de Sines, empresas que nele 
participam, planeamento das instalações e investimentos pre- 
vistos. 

Descreve-se depois, de uma forma sumária, O estado dos 
trabalhos em curso, em fins de Agosto de 1977, 
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U, D. C. 621,39(469) 


MANUEL JOSE ABREU FARO 
PEDRO LUIS BORGES TEIXEIRA 


CONCEPTION AND PLAN OF TELECOMMUNICATION 
SYSTEMS IN PORTUGAL 


Técnica N.º 445/446 — LIV LH] — 12/1-1977/78, p. 275 a 295 


U. D. C, 665.51/58 
LUIS L. RAMALHETE 
PETROCHEMICALS 


Técnica N.' 445446 — LIWLII — 1204977/78, p. 297 a 303 


U. D. C. 665.51/.58 
JOAO PEDRO ANTUNES GOMES 


THE SINES PETROCHEMICAL COMPLEX; THE CONCEPTION 
AND ACTUAL DEVELOPMENT OF THE WORKS 


Técnica N.º 445/446 — LII/LII — 12/1-1977078, p. 305 a 312 


After mentioning the legal documents that created the 
olefins petrochemical industry, we describe the conception of 
the Sines Petrochemical Complex, participating companies. 
installations planning and foreseen investments. 

Afterwards, we summarise the actual development of the 
works, in the end of August 1977, 


U. D, C. 624.07[012.4]8 
S. POMFEU DOS SANTOS 
CONNECTIONS IN PRECAST CONCRETE STRUCTURES 
Técnica Nº 445446 — LILI — 121407778, p. 259 a 254 


Precast building structures demand that connections bet- 
ween precast elements be made on the job, Design, project 
and production of connections require the resolution of certain 
specific problems, which are not usually envisaged in regula- 
tions for reinforced end prestressed concrete structures and 
steel structures. 

In this paper connections between precast elements of 
concrete framed structures are dealt with from the point of 
view of such requirements as structural safety, durability 
and handling, in order to show how these requirements may 
be fulfilled. A classification of connections according to some 
parameters |s also presented and some considerations ere 
mede about design and execution of some types IF connections. 
Several examples of connections made by different methods 
are given. Finally the author presents a programme of studies 
to be undertaken within this scope. 


U, D. C. 627.8.09(679.4) 


ALBERTO ABECASIS MANZANARES 
ANTONIO CARRISSO 

ANTÔNIO DE CARVALHO QUINTELA 

ANTÔNIO JOSÉ DA COSTA SANTOS 
JOÃO DE SALVADOR FERNANDES 


PORTUGUESE ENGINEERING AND THE INSTALLATION 
OF LARGE POWER UNITS, THE CASE OF CABORA-BASSA 


Técnica N.º 445/4468 — LIL/LII — 1211-1977/78, p. 255 a 266 


The hydro-electric scheme of Cabora-Bassa is presented 
in order to attest the capacity of the Portuguese Engineering 
in the concsption and design of power stations equipped with 
large units. 

Afier a general description of lhe scheme, the principal 
constrainis and problems which have been confronted are 
brigfly referred, in order to emphasise the necessity that they 
must be joinily solved by the equipment manufecturers and 
the Consultents. 


U. D. C. 621,313.12:519.,28 


M. FERREIRA DE OLIVEIRA 
1. REIS CARDOSO 


SPINNING RESERVE STUDIES INCLUDING GENERATOR 
OPERATED STUDIES 


Técnica N.' 445/446 — LIWLIHH — 121497778, p. 267 a 274 

Spinning reserve studies including reliability techniques 
allow the risk evaluation, i.e. the probability of the system 
load exceeding the system capacity. 

This studies require the knowledge of the individual pro- 
babilities of each of the generator states; catch state is de- 
fined by the maximun capacity the generator is able to pro- 
duce. 

This paper presents a technique to evaluate each of the 
individual probabilíties, as a function of time, of the generator 
states required in spinning reserve studies. 
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C. D. U. 622.792+ 622.798 
H. MAIA E COSTA 


AGLOMERAÇÃO DE CONCENTRADOS FINOS 
DE MINÉRIOS 


Técnica N.º 445/4465 — LIILIHI — 121-1977/78, p. 349 a 354 


Nos modernos processos da metalurgia extractiva exigem- 
-se cargas devidamente apropriadas de modo a conseguirem-se 
produtividades tão elevadas quanto possivel e custos ener- 
géticos tão reduzidos quanto possivel. 

Neste trabalho indicam-se os processos de aglomeração 
mais difundidos mundialmente e que tem a sua maior apli- 
cação na indústria siderúrgica: sinterização e peletização. 

Trata-se depois de inventariar o que já existe em Portugal 
e quais as perspectivas futuras tendo em atenção os recursos 
nacionais em minérios de ferro. 


C. D. U. 669.1:543,42.64 
M. HELENA DE MELO CARVALHO 


DETERMINAÇÃO DE COEFICIENTES DE REPARTIÇÃO 
SOLIDO/LIQUIDO EM LIGAS TERNÁRIAS Fn-Mn-C 


Técnica N.º 445/446 — LILI — 121 4977/78, p. 3255 a 364 


Usando a técnica de análise por microsonda, determi- 
naram-se coeficientes de repartição sólido/liguido para algu- 
mas ligas da região rica em Fe do sistema Fe-Mn-C. O mé- 
todo utilizado envolveu a microanálise de amostras tempe- 
radas a partir de temperaturas de equilibrio na região L+%. 
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C. D. U. 621.3.016.33:518.5 


1. M. FERREIRA DE JESUS 
J. P. SUCENA PAIVA 


SOLUÇÃO DO TRÂNSITO DE POTÊNCIA PELO METODO 
DO DESACOPLAMENTO 


Técnica N.º 445/4465 — LI|/LIHI — 121-4977/78, p. 355 a 382 


O cálculo do trânsito (ou fluxo) de potência numa rede 
de energia eléctrica, é presentemente efectuado em computa- 
dor digital por meio de diversos métodos numéricos, dentre 
os quais o mais sofisticado é o método de Newton-Raphson 
e suas variantes. 

No artigo os autores descrevem um programa de cálculo 
baseado numa destas variantes, o método do desacoplamento, 
inicialmente proposto por Stott [10]. 

E dada especial ênfase és técnicas de programação com 
esparsidade utilizadas. Apresentam-se os resultados numéricos 
obtidos para diversas redes, comparando-as com os obtidos 
pelo método de Newton-Raphson clássico. 


C. D. U. 622.2:519.3 
JOSE QUINTINO ROGADO 


OPTIMIZAÇÃO DO BENEFICIAMENTO DE SISTEMAS 
DE JAZIDAS MINERAIS. APROPRIAÇÃO RACIONAL 
DOS RECURSOS MINERAIS 


Técnica N.º 445/4466 — LILI] — 12/1-1977/78, p. 383 a 388 


No planeamento a níveis nacional ou regional e também 
no âmbito de uma grande empresa responsável por vários 
depósitos minerais, é-se conduzido ao problema da optimiza- 
ção global de sistemas de jazidas, os quais podem conter 
minérios diferentes. 

Tomando como constrangimentos a taxa anual de produ- 
ção de minério tal-qual e o teor minimo dos concentrados 
aceite pelos utilizadores, o objectivo consiste em optimizar 
duas funcionais simultaneamente: a perda total em metal e a 
margem monetária global não actualizada. 

Para rresolver este problema, usa-se a «função caracte- 
risticas, já apresentada noutra oportunidade, mediante um 
processo de «homogenizaçãos que se traduz na utilização de 
«teores reduzidos». 

Deste modo, constroi-se a função caracteristica global 
do sistema e pode obter-se facilmente a margem monetária 
optimizante que conduz à perda minima em metal compativel 
com os constrangimentos. 

O programa CARACT que optimiza as duas mencionadas 
funcionais é sumarizado e algumas das suas virtualidades 
analisadas. 
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U. D. C. 621.3.016.33:518.5 


J. M. FERREIRA DE JESUS 
J. P. SUCENA PAIVA 


POWER FLOW SOLUTION BY THE DECOUPLING METHOD 
Técnica N.º 445/446 — LILI — 12/1-1977/78, p. 365 a 382 


The power flow solution in an electric eneroy system is 
currently executed by digital computer using several numerical 
methods, of which the most powerful is the Newton-Raphson 
and its variations, 

In this paper the authors describe a computer programme 
based upon one of these variations, the decoupling method, 
initially suggested by Stott [10]. 

The sparsity programming thecniques are specially empha- 
zised. Numerical results for several real systems are presen- 
ted, and compared with those obtained by the classic Newton- 
-Raphson method. 


U, D.C. 622,2:519.3 
JOSE QUINTINO ROGADO 


RATIONAL APPROPRIATION OF MINERAL RESQURCES AN 
OPTIMIZING CRITERION FOR BENEFICIATION OF SYSTEMS 
OF ORE DEPOSITS 


Técnica N.º 445/446 — LII/LII — 12/1-4977/78, p. 383 a 388 


When planning for large areas, say a country, à mining 
district, or a great concern responsible for several mines, 
the exploitation and beneficiation polícies deal with complex 
avetems, 

in à global view point, the best objective |s tô minimize 
the overall loss of metals and to maximize the total non- 
-discounted monetary margin. 

The most important constraints framing the resolution of 
such a problem are the minimum (or the maximum) grade 
levels of the concentrates to be produced and the annual 
rates of production admissible by the market. 

The method developed in this paper is a quite general 
one and uses the concept of characteristic function elsewhere 
developed by means of an homogenization process which leads 
to consider the deposits set as a sole ore body, even If far 
apart ore bodies there exist. 

The computer program CARACT (built up for the am 
method) is outlined and results already obtained with it are 
summarized. 

Finally, among other suggestions, the possibility of using 
this method for the evaluation of mineral resources In world 
scale, |s advanced, 


U. D. C. 622.792 + 622.788 
H. MAIA E COSTA 
AGGLOMERATIOMN OF FINE ORE CONCENTRATES 
Técnica N.º 445/4466 — LI|LII| — 1214077/78. p. 340 a 354 


in extractive metallurgy modern processes for preparing 
tha charges are required in order to reach as high as possible 
productivity and as low as possible energetic costs. 

In this work are indicated the more worldly widspread 
agglomeration systems and those that have their main appli- 
cation in the siderurgical industry: sintering and pelletizing. 

After that an inventory of what already exists in Portugal 
is drawn up together with its future prospects, paying atten- 
tion to the national resources in iron ore. 


U. D. €. 669.1:543.42.64 
M. HELENA DE MELO CARVALHO 


DETERMINATION OF SOLIDLIQUID PARTION 
COEFFICIENTS IN Fe-Mn-C TERNARY ALLOYS 


Técnica N.º 445/446 — LILI] — 12H A97TI7B. p. 355 a 364 


Solidiliquid partition coefficients of some Fe-Mn-C alloys, 
at the Fe — rich corner, have been determined by means of 
electron probe microanalysis. The method used involved the 
microanalysis of specimens quenched from an equilibrium 
temperature In the L+4y field. 
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C. D. U. 669.14:621.785:620.163 
ANTERA VALERIANA DE SEABRA 


O CONTROLE DE QUALIDADE NOS TRATAMENTOS 
TÉRMICOS 


Técnica N.º 445/446 — LII/LHI — 12/1-1977/78, p. 313 a 326 


Em face da importância cada vez maior do «Controle de 
Qualidade» para garantia do nivel de qualidade das peças 
obtidas segundo dado processo de fabricação referem-se no 
|| Capitulo algumas considerações fundamentais, relativamente 
a peças sujeitas a tratamento térmico desde a selecção da 
matéria-prima até ao controle do produto final. 

Pára um melhor enquadramento do assunto recordam-se 
num primeiro capitulo as noções gerais relativas aos conceitos 
de Qualidade, Controle e Controle de Qualidade. 


C. D. U. 621.762:.689.1 
CARLOS ALBERTO PACHECO DA SILVA 
PULVEROMETALURGIA DE LIGAS FERROSAS 
Técnica N.º 445/4465 — LIF/LIHN — 12/1-1977/78, p. 327 a 338 


Após se realçar o interesse da pulverometalurgia como 
tecnologia alternativa de enformação, faz-se uma breve des- 
crição das etapas do processo e o modo como podem in- 
fluenciar as propriedades. 

Referem-se finalmente os tipos de ligas ferrosas mais 
frequentemente produzidas por pulverometalurgia, 
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Ano Lil | LIIHI — Dezembro 1977 | Janeiro 1978 


C. D. U. 621.791:621.311.2 
SALGADO PRATA 


ALGUNS ASPECTOS DA TECNOLOGIA DE SOLDADURA 
NA CONSTRUÇÃO DE EQUIPAMENTOS PARA PRODUÇÃO 
DE ENERGIA 


Técnica N.º 445/4466 — LII/LIH — 121497776, p. 339 a 343 


Trata o tema Construção no que respeita à tecnologia de 
soldadura, aplicada a equipamentos para centrais térmicas, 
hidráulicas e nucleares, tal como actualmente se apresenta 
nã empresa em que o autor exerce a sua actividade. 

Faz-se um breve inventário dos processos de soldadura 
utilizados pela empresa na construção de equipamentos de 
energia. 

Indicam-se limitações dos processos de soldadura pro- 
venientes do risco de provocarem deterioração durante o 
fabrico (fissuração a quente, a frio, durante tratamentos tér- 
micos) ou em serviço, por exemplo, ruptura frágil ou corrosão. 

Dascrevem-se e exemplificam-se algumas técnicas para 
rodear as limitações mencionadas. 

Referem-se alguns ensaios em cujos resultados se baseiam 
as técnicas adoptadas, 


C. D. U. 621.785.001.2:869.14 
ANTÔNIO DE PADUA LOUREIRO 


NECESSIDADE DE INVESTIGAÇÃO E DE DESENVOLVIMENTO 
TECNOLÓGICO EM PORTUGAL NA AREA DOS TRATAMEN- 
TOS TÉRMICOS DOS AÇOS 


Técnica N.º 445446 — LILI] — 12/1-1977/78, p. 345 a 348 


Portugal importa & quase totalidade dos aços de trata- 
mento térmico que emprega na indústria. À falta, na maioria 
das empresas, de instalações de tratamentos térmicos devida- 
mente apetrechadas e de pessoal com os indispensáveis 
conhecimentos, acarreta a rejeição de muitas peças logo após 
o tratamento ou a sua curta duração em serviço, causando 
grandes perdas para o Pais em divisas e em mão-de-obra. 
Por outro lado, a falta de investigação suficiente para reso- 
lução dos problemas que a indústria portuguesa defronta, não 
permite a racionalização dos tipos de aços a importar, impe- 
dindo uma decisão fundamentada sobre a eventual criação 
de uma acearia de aços especiais em Portugal, nem permite 
a optimização dos tratamentos mais adequados a cada caso. 

Na presente comunicação apontam-se os problemas mais 
importantes nesta área, sugerindo-se seguidamente a criação 
de uma associação técnica e de uma oficina nacional de 
tratamentos térmicos dos aços, com vista a incentivar uma 
progressiva independência tecnológica de Portugal neste cam- 
po e a suprir as deficiências em instalações e em pessoal 
especializado com que lutam muitas empresas que necessitam 
de efectuar tratamentos térmicos a peças que produzem. 
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U, D. C. 621.791:621.311.2 
SALGADO PRATA 


SOME FEATURES OF WELDING TECHNOLOGY IN THE 
POWER EQUIPMENT FABRICATION 


Técnica N,” 445446 — LI|/LIHI — 121497778, p. 38 a 38 


lt deals with the subject construction on what is con- 
cerned with welding technology, applied to equipments for 
thermal, hydraulic and nuclear power plants as now it is 
presented in the department where the author works. 

A brief inventory of the welding processes used by the 
factory ou the making of energy equipments is done. 

Limitations of the welding processes because of the danger 
of causing deterioration during the making (fissures during 
thermal treatments) or on usage, for instance, fragile rupture 
or corrosion are indicated. 

Some technics to dodge the above limitations are des- 
cribed and examples are given. 

Reference is made to some experiments on which results 
the adopted thecnics are based. 


U. D. C. 621.785:001.2:660.14 
ANTÔNIO DE PADUA LOUREIRO 


HEAT TREATMENT OF STEELS: RESEARCH AND 
DEVELOPMENT NEEDS IN PORTUGAL 
Técnica Nº 445/446 — LIVLIH — 121-41977/78, p. 345 a 348 

Portugal imports nearly all the seat treatablo steels needed 
by her industry. The great majority of our factores lack 
suitable facilities and technical staff with appropriate expe- 
rence on heat treating steel, This situation leads to the 
rejection of many parts which were just heat treated or upon 
a short period under service conditions. Spending in foreign 
currency end labouyr costs is thus unavoidable. On the other 
hand, insufficient research on industrial problems elther thus 
not allow a sensible selection of the steel types which should 
be imported thus affecting a well based decision on the 
creation of a special steels industry or reduces the ability 
to optimization In each case. 

In this paper a review is made of the most important 
problems in this area. Two suggestions are put forward: the 
creation of a technical association and the sitting up of à 
steel heat treatment work shop at a national level. In both 
cases the end results almed at are: 1) to make our country 
progressively less technologically dependent on foreign know 
how and imports: end 2) to give technical support to those 
firms which have neither appropriate facilities nor sufficiently 
tralned personnel to deal with their heat treatment problems. 


U, D. C. 660,14:621.785:620.163 
ANTERA VALERIANA DE SEABRA 
QUALITY CONTROL IN HEAT-TREATMENT 
Técnica N.º 445/4468 — LII/LIN — 12/1-1977/78, p. 313 a 326 


Taking into account the increasing importance of Quality 
Control to guarantee the level of quality of parts obtained by 
any manufacturing process, in Chapter |l some basic conside- 
rations are made as regards heat-treated parts, starting from 
the selection of raw-materials up to the control of the final 
product. 

To put the matter inte its right perspective, Chapter | 
recalle some general notiens on the concepts of Quality 
Control and Control of Quality. 


U. D. C. 621.762:669.1 
CARLOS ALBERTO PACHECO DA SILVA 
IRON POWDER METALLURGY 
Técnica N.º 445/4468 — LII/LIN — 12/-1977/78, p. 327 a 338 


After emphasizing the interest of powder metallurgy as 
an alternative technology of metal forming, the steps of the 
process are described together with the way they influence 
the properties, 

Finally the more usual ferrous alloys produced by powder 
metallurgy are referred. 
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“À nossa vocação é exportar 


Caixa de escape para uma central nuclear, 
destinada à Electricité de France 


EQUIPAMENTOS HIDROMECÂNICOS 


EQUIPAMENTOS ELECTROMECÂNICOS 
PARA CENTRAIS HIDROELECTRICAS 

E TERMICAS 

(CLÁSSICAS E NUCLEARES) 


EQUIPAMENTOS PARA A INDÚSTRIA 
QUÍMICA E PETROLÍFERA 


MATERIAL CIRCULANTE DE CAMINHO 


DE FERRO Países para onde temos exportado : 
EQUIPAMENTO OFFSHORE PARA AFRICA DO SUL - ANGOLA - ARGENTINA - BRASIL - COSTA 
PROSPECÇÃO E EXPLORAÇÃO RICA - FRANÇA - GUATEMALA - IRÃO - IRAQUE - ISLÂNDIA 
DE PETROLEO KUWEIT - MALAWI : MARROCOS - MÉXICO - MOÇAMBIQUE 
CONSTRUÇÕES MODULADAS, NICARÁGUA - PAQUISTÃO - RODÉSIA - SUDÃO : TURQUIA 
PONTES E ESTRUTURAS UGANDA - U.S.A. - VENEZUELA - ZAIRE 


SOREFAME 


SOCIEDADES REUNIDAS DE FABRICAÇÕES METÁLICAS, S.A.R.L. 
Amadora — Portugal 


À procura da água (mole) 
ou da pedra (dura) 


À procura de qualquer formação rochosa, independentemente 
da sua consistência — mole ou dura — a linha ROC 
actua com a maior eficiência. 


Seja onde for o seu trabalho de captação de água, escavação 
terraplanagem, sondagem, ancoragem, abertura de túneis, etc,, 
a Atlas Copco possui o carro de perfuração adequado. 
Manobrável por um só homem. 


De facto, a linha ROC, a mais vasta gama de carros 

de perfuração sobre lagartas, equipados com colectores de 
poeiras, permite-lhe abrir furos de 35 a 216 mm de diâmetro 
a profundidades para além de 200 m, 

com um rendimento insuperavel. 

No entanto, se o seu caso exige outro equipamento 

de perfuração, de capacidade inferior ou superior, 

a solução está ao seu alcance. Exponha-nos o seu problema. 


Contacte o Centro de Ar Comprimido mais próximo. 
Soc. Atlas Copco de Portugal, Lda. 
Centros de Ar Comprimido: 


Lisboa (Linda-a-Velha) — Av. 25 de Abril — Tel. 2190011 
Porto — Rua Eng.º Ferreira Dias, (à via rápida) - Tel. 691071 
Vila Viçosa — Largo D. João IV, 60 —- Tel. 42101 


